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ABSTRACT 

We present the far-IR, millimeter, and radio photometry as well as optical and near-IR spec- 
troscopy of a sample of 48 Spitzer-selected ULIRGs with IRS mid-IR spectra. Our goals 
are to compute their bolometric emission, and to determine both the presence and relative strength 
of their AGN and starburst components. We find that strong-PAH sources tend to have higher 
MIPS160/im and MAMBO 1.2 mm fluxes than weak-PAH sources. The depth of the 9.7 /xm sihcate 
feature does not affect MAMBO detectability. We fit the far-IR SEDs of our sample and find an 
average (im) 7 x 10^^ Lq for our z > 1.5 sources. Our spectral decomposition suggests that strong- 
PAH sources typically have ~20-30% AGN fractions of Lir . The weak-PAH sources by contrast 
tend to have ^ 70% AGN fractions, with a few sources having comparable contributions of AGN and 
starbursts. The optical line diagnostics support the presence of AGN in the bulk of the weak-PAH 
sources. With one exception, our sources are narrow-line sources, show no obvious correspondence 
between the available optical extinction and the silicate feature depth, and, in two cases, show some 
evidence for outflows. Radio AGN are present in both strong-PAH and weak-PAH sources. This is 
supported by our sample's far-IR-to-radio ratios (q) being consistently below the average value of 
2.34 for local star-forming galaxies. We use survival analysis to include the lower-limits given by the 
radio-undetected sources, arriving at (g) = 2. 07 ±0.01 for our z>1.5 sample. In total, radio and, 
where available, optical fine diagnostics support the presence of AGN in 57% of the z > 1.5 sources, 
independent of IR-based diagnostics. For higher-z sources, the AGN luminosities alone are estimated 
to be >1O^^L0, which, supported by the available [Oiii] luminosities, implies that the bulk of our 
sources host obscured quasars. 

Subject headings: galaxies:infrared, radio, optical AGN 



1. INTRODUCTION 

Ultraluminous Infrared Galaxies (ULIRGs; defined as 
having Lir,>1O^^L0, (Soifer et al. 1987; Sanders & 
Mhabel 1996; Genzel et al. 1998; Lonsdale et al. 2006)) 
are believed to be an important contributor to the global 
star-formation rate density at z^2 (Le Floc'h et al. 
2005; Caputi et al. 2007). This simple luminosity-cut 
definition, however, potentially includes a wide range 
of physical characteristics. It is particularly poorly un- 
derstood at the high luminosity (Lir ~ 10^^'^ - 10^^ L©) 
end, due to the scarcity of such sources locally. In this 
regime and above, most sources are believed to be largely 
AGN-dominated (Lutz et al. 1998; Tran et al. 2001; Far- 
rah et al. 2002). However, at z^2, sub- mm galaxies 
(SMGs) fall in the above luminosity range, but appear 
to be largely starburst-dominated ULIRGs (Alexander 
et al. 2005; Vahante et al. 2007; Pope et al. 2008). Al- 
though, less numerous, luminous quasars also peak at 
z ~ 2 (Richards et al. 2006) , and their bolometric lumi- 
nosities fall in the above luminosity range. A popular 
evolutionary scenario suggests that major mergers lead 
to starburst-like cold ULIRGs followed by warm ULIRGs 
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as the dust-obscured AGN become dominant. Eventu- 
ally, outflows shed the dusty cocoon to reveal the naked 
quasar (Sanders et al. 1988). Due to their mid-IR bright- 
ness Spitzer Space Telescope (Werner et al. 2004) se- 
lected z~2 ULIRGs, are biased to the 'warm ULIRG' 
phase (Yan et al. 2004b, 2005; Houck et al. 2005) and so 
fill the gap between SMGs and QSOs allowing us to test 
the above scenario directly at the peak epoch for ULIRGs 
and quasars. If at least a part of the quasars follow a sim- 
ilar evolutionary path we expect to find a population of 
obscured quasars which are accompanied by strong star- 
burst activity. Recent studies suggest this might indeed 
be the case (Polletta et al. 2008; Martfnez-Sansigre et al. 
2008). 

This paper concerns a sample of 52 Spitzer-selected 
ULIRGs, with Infrared Spectrograph (IRS; Houck et 
al. 2004) spectra. By selection, our sample has some- 
what less extreme 24 fj,m -to-R colors compared with 
other samples aiming at understanding the Spitzer- 
bright (> 1 mjy at 24 /^m), optically- faint population 
(e.g. Houck et al. 2005; Weedman et al. 2006). There- 
fore, beside the strong continuum, strong silicate absorp- 
tion spectra typical of such samples, our sample also in- 
cludes sources with PAH features of comparable strength 
to those seen in SMGs (Pope et al. 2008). Our sample 
may therefore represent transition objects spanning the 
range between (un)obscured quasars and SMGs. This 
is Paperlll of a series of papers on the infrared proper- 
ties this sample. In Yan et al. (2007, hereafter Paper I), 
we presented the spectra, sample selection as well as the 
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IRS-derivcd rcdshifts. In Sajina et al. (2007a, hereafter 
Paperll), we presented the spectral decomposition of the 
IRS spectra into their PAH, continuum, and obscuration 
components. 

An HST NICMOS morphology study of the 36 ^ > 1.5 
sources in our sample is presented in Dasyra et al. (2008). 
That study shows that, typically, the most compact 
sources are those with AGN-like mid-IR spectra (i.e. 
weak PAH) and vice versa. The more intriguing re- 
sult was the finding that at least some of the interac- 
tions leading to ultraluminous activity at 2;~2 appear 
to differ from local ones (which essentially require equal 
mass mergers). This result cautions that some of our 
sources may not have true local analogues. In Sajina et 
al. (2007b), we discussed the fact that 40% of our z>1.6 
sources with optical depths of Tg 7 > 1 are radio loud, in- 
cluding some sources with powerful radio jets. Since deep 
silicate feature sources are particularly prone to power 
source ambiguities (due to the heavy obscuration), this 
result was crucial in confirming that a significant fraction 
of these sources do indeed host powerful AGN. Lastly, the 
Chandra X-ray detcctability of our sources is addressed 
in Bauer et al. (2008, in prep.). There are however 
key questions about the nature of these sources which 
remain unanswered after the above studies, and are the 
subject of this paper. Given the uncertainties inherent in 
the mid-IR diagnostics (see Papcrll) , can wc confirm the 
presence of obscured AGN in our sources, and if present, 
are they of quasar-like strength? What is the bolometric 
power output of our sources, and what fraction of that 
is due to AGN vs. star- formation activity? 

This paper is organized as follows. § 2 presents the 
available multiwavelength data and data reduction in- 
cluding Keck and Gemini optical and ncar-IR spec- 
troscopy, Multiband Imaging Photometer for Spitzer 
(MIPS; Rieke et al. 2004) 70 /xm and 160 /itm, and 
MAMBO 1.2 mm pointed photometry, as well as the 
archival VLA 1.4 GHz, and GMRT 610 MHz data. Our 
analysis of the above data is presented in § 3, where we 
look for the signature of AGN activity via their radio 
properties (luminosity and spectral index) as well as tra- 
ditional optical line diagnostics (broad lines, high ion- 
ization lines, and AGN-like line ratios). We use spec- 
tral decomposition to determine the role of AGN and 
starburst activity in powering the bolometric output of 
our sources. In §4, we use all of the above diagnos- 
tics to address the qiiestion how do our sources relate to 
other high-z luminous populations such as SMGs and in 
particular whether or not they form part of the popu- 
lation of obscured quasars as we speculated in Paperl. 
Finally, we summarize the key results of this work in 
§ 5. Throughout this paper, we adopt a ACDM model 
with f]M=0.27, f}A=0.73, and ffo = 71 kms'^Mpc"^ 
(Spergel et al. 2003). 

2. DATA 

2.1. IRS spectral sample 

Paperl gives the detailed description of the sample se- 
lection, presents the observed spectra and redshift mea- 
surements. Here we summarize the key points. Our 
sample is selected from the 4 sq.deg. Spitzer Extragalac- 
tic First Look Survey (XFLS*^). The full sample con- 
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sists of 52 sources with 24 /im flux density brighter than 
0.9mJy, and red 24-to-8/im and 24-to-R colors^. The 
IRS low resolution data revealed that 47 of the targets 
have measurable redshifts from the features detected in 
their mid-IR spectra, with the majority (35/47 = 74%) 
at 1.5 <z<3 .2, and a smaller fraction (12/47 = 28%) at 
0.65<z<1.5. An additional source (MIPS279) has a 
Keck redshift bringing our sample to 48. Most of our 
spectra cover rest-frame ~ 5 - 15 iim with a minimum 
signal-to-noise per pixel of 4. Because the primary se- 
lection is the 24 /xm flux, we name our sources based on 
their MIPS24/xm catalogue number (see Table 1 for a 
source list and redshifts). 

In Paperll, wc presented the decomposition of the 
rest-frame 5- 15 fim mid-IR spectra into PAH, contin- 
uum, and obscuration components. Based on the 7.7 /um 
PAH feature equivalent widths, we define 'strong-PAH' 
sources as those having EW7.7 > 0.8/im , while 'weak- 
PAH' sources have EW7.7<0.8/xm. Wc adopt these 
definitions throughout this paper as well. In Paperll, we 
found that ~ 75% of the sample are weak-PAH sources 
consistent with AGN-dominancc in the mid-IR. How- 
ever, although of comparable luminosities, these sources 
are more obscured than typical quasars as evidenced 
by their steeper slopes, and typically significant 9.7 /um 
silicate absorption. Although most of our strong-PAH 
(starburst-likc) sources were found among the lowcr-z 
sub-sample, we found an intriguing population of these 
star-forming sources at 2 ~ 2 with luminosities an order 
of magnitude larger than local ULIRG analogues. 

2.2. Optical spectroscopic observations 

We obtained optical spectra for 12 of our sources using 
the LRIS and DEIMOS instruments on the Keck tele- 
scopes (Oke et al. 1995; Faber et al. 2003). The initial 
results from the optical and near-IR spectroscopy (see 
below) were presented in Paperl (Yan et al. 2007). For 
the red side of the LRIS spectrograph, he D680 dichroic 
was used, splitting the red and blue arms at 68OOA. The 
red arm was equipped with the 400 line mm~^ grating 
blazed at 8500-lOOOOA. At the blue side, we used the 
300 line mm^^ grating blazed at 5000A. Roughly, the 
full spectral range is 3200-9600A. The slit width is 1.2 
and we dithered the targets along the slits between ex- 
posures with offsets of 2-4 . The integration time for 
each exposure was 1200 seconds, with a total time per 
mask of 1-1.5 hours. The DEIMOS observations were 
made using a 1200 line mm~^ grating with central wave- 
length settings of 7400A and 7699A together with the 
GG495 blocking filter, resulting in a O.33A pix ^ mean 
spectral dispersion and a 1.45A instrumental resolution. 
The wavelength coverage for DEIMOS spectra is 6300- 
93OOA. The total integration time for DEIMOS observa- 
tions is 3,600 seconds per mask. 

For the LRIS data reduction, we used a suite of IRAF 
scripts written by D. Stern (private commmiications). 
The basic procedure is described in detail in Yan et al. 
(2004a). For the DEIMOS data, we used the DEEP2 
SPEC2d data reduction pipeline, which performs cosmic 
ray removal, flat-fielding, co-addition, sky subtraction, 

(1) R(24,8) = logio(j-/,.(24/im)M,(8/im) > 0.5; (2) 
i?(24,0.7) = logio{uU{24ij,m)/uf^{0.7ij.m) > 1.0. 
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wavelength calibration and both 2-d and 1-d spectral ex- 
traction. 

2.3. Near-IR spectroscopic observations 

We obtained ncar-IR H and K spectra using the NIR- 
SPEC instrument on Keck (McLean et al. 1995). Addi- 
tional three nights with the NIRI instrument (Hodapp 
et al. 2003) on the GEMINI-north telescope were used 
to obtain near-IR H and -Pf-band spectra of our sample 
targets at z ^ 1.8 — 2.4. Ncar-IR spectra with suffi- 
cient signal-to-noise ratios were obtained for 11 sources. 
The NIRSPEC observations were done in moderate res- 
olution of ^ 1500, with a slit of 42 in length and 
0.76 (4pixels) in width. At this resolution, the NIR- 
SPEC instrument FWHM is ~ 10.0 A(~ 200 km/s) at H 
and ^ 14A(~ 190km/s), as determined from isolated sky 
lines (sec also Brand et al. 2007). The NIRI observa- 
tions wore taken with the f/6 camera, with a field of 

view of 120" X 120". We selected the slit width of 0.75" 
(6.42 pixels) and achieved roughly 520 spectral resolu- 
tion in H (1.43-1.96 /im) and K-band (1.90- 2.49 /im). 
The NIRI instrument FWHM is ~ 20 A(~ 400km/s) at 
H and ~ 28 A(~ 380km/s) in K (Brand et al. 2007). The 
soiirccs for the ncar-IR spectroscopy were chosen to have 
appropriate redshifts so that their expected Ha and [Oiii] 
emission lines are within the observed H or ii'-band. The 
total on-target integration time per source ranges from 
45minutes to 1.5 hours, depending on the strength of the 
emission lines and the observing conditions. 

2.4. Near-IR data reduction and flux calibration 

The NIRSPEC data reduction was done using IDL 
software written by George Becker (private communi- 
cation). This software uses the methods of optimal sky 
subtraction (Kelson 2003) for reducing longslit spectra in 
the near- infrared, where the strong OH emission lines can 
hide the faint spectral features of our source sample. This 
method models and removes sky emission without first 
rectifying or rcbinning the data, which eliminates noise 
and edge effects that can occur during interpolation. The 
spectra are wavelength calibrated using the sky emission 
lines. Standard F and G stars were observed in order to 
perform flux calibrations and to remove atmospheric ab- 
sorption features. The absolute standard fluxes were ob- 
tained from the IRTF Spectral Library^ using the spec- 
tral type closest to that of our standards. The maximum 
air mass difference between our source observations and 
the standard star observations is ~ 0.2. Given the av- 
erage values of sky extinction on Mauna Kea, which are 
0.05mag/airmass, and 0.07mag/airmass in the H and K 
band respectively (Krisciunas et al. 1987), this airmass 
difference does not significantly affect our flux calibra- 
tion. 

The aforementioned IDL code for longslit reductions 
does not currently accept NIRI data, so we performed 
the spectroscopic reductions on the NIRI data using 
the Gemini IRAF package GNIRS and the instrument- 
specific NIRI package, following the data reduction pro- 
cedures outlined on the Gemini Observatory webpages^. 

® http:/ /irtfwcb.ifa. hawaii.edu/ spcx/spcxlibrary/IRTFlibrary.html 
^ http;//www.gemim.edu/sciops/instruments/niri/data-format- 
and-reduction 



The flux calibration followed a similar procedure as for 
NIRSPEC. The spectra were extracted using a Gaussian 
profile with FWHM of 4pixels for NIRSPEC and 8pixels 
for NIRI, appropriate for spatially unresolved sources. 

Lastly, the absolute flux correction is uncertain due 
to the uncertain positioning of the standard sources on 
the slits and possible light loss (see also Brand et al. 

2007) . In order to estimate the degree of uncertainty 
in the absolute fluxes, we compare with expectations 
from the NICMOS broadband photometry (Dasyra et al. 

2008) . This is complicated for most of our sources due 
to the non-detection of the continuum. The only excep- 
tion is MIPS8196 (observed with NIRI), which is brighter 
than the rest of our sources {mH,Vcga = 18.68). We find 
that, convolving its spectrum with the NICMOS F160W 
transmission curve results in a flux 1.6 times the ob- 
served. This sources also has a K -hand magnitude from 
the Wide Field IR Camera (WIRC), of niKyegi, = 17.63 
(P.Choi, private communication). Our if -band spectrum 
results in a WIRC fltix a factor of 1.9 higher than the ob- 
served. This shows that the absolute flux calibration of 
our sources is likely overestimated by roughly 50 100%. 
However, as we cannot confidently determine these over- 
estimation factors for any other sources (both NIRSPEC 
and NIRI), we choose not to apply this correction factors, 
but merely state that the line luminosities we derive from 
these spectra (§ 3.4.2 and Table 4) are uncertain up to a 
factor of 2 (0.3 dcx). 

2.5. MIPS lOiim 

The XFLS was observed in scan map mode with the 
MIPS70/im array (Fraycr et al. 2006a). A smaller, 0.25 
sq.deg., verification field was imaged with 4 times the in- 
tegration time as the main XFLS field. The noise is non- 
uniform across the field, however the average point source 
(Icr) rms at 70 /xm is ~2.8mJy in the main field, and 
~ 1.6 mjy in the verification field. Three of our sources 
fall inside the verification field (MIPS15977, MIPS8268, 
MIPS8245). For the sources without clear 70 /im detec- 
tion in the scan map data, additional pointed photome- 
try data was obtained as part of Spitzer G02 program 
#20239 (PI: P.Choi). Most of our targets have roughly 
700 seconds integration time per pixel. The scan map 
data and pointed photometry data were combined into a 
local median-subtracted image with 4" pixel scale. 

The fluxes and their uncertainties are estimated from 
PRE photometry on the median-subtracted image (using 
APEX; Makovoz & Marleau 2005). For color corrections, 
we use the values for flat spectrum {F„(xi'~^), which 
is also roughly consistent with the color corrections for 
dust temperatures in the range 30-70K. This implies 
a color correction of 1.09 for MIPS70/im, and 1.04 for 
MIPS160/zm (see below). An additional 5% error in the 
flux is attributed to the uncertain color correction (after 
excluding obviously inconsistent spectral types). This 
uncertainty is negligible given the generally low signal-to- 
noise of our data. The typical rms flux error is ^ 1.7 mJy. 

We detected 29 of the 48 sources listed in Table 1, 
based on a 10" matching radius (the positional off- 
sets and possibility of random matches are addressed in 
§2.11). The aperture-photometry derived 70 /zm fluxes 
for our sample were presented in Paperll. We found that, 
within the expected uncertainty the PRF-derived fluxes 
presented here and these aperture-photometry derived 



4 



Sajina et al. 



fluxes agree well, however, that the PRF method is some- 
what more conservative in that 4 sources with low signal- 
to-noise detection in Paperll (SNR'^S-S.S) were not de- 
tected with APEX. Here we adopt the PRF fluxes as 
being more conservative. 

2.6. MIPS 160 ^im 

The XFLS field and verification field were also im- 
aged with the MIPS160/im array (Frayer et al. 2006a). 
The noise is quite non-uniform but typically 20 mJy 
in the main field, and ~ lOmJy in the verification field. 
Two of our sources are in the MIPS160/im catalogue: 
MIPS22404 and MIPS15977 (which is in the verification 
field). Our sample, except MIPS224G4, was also observed 
in 160 /im pointed small-field photometry {Spitzer G02 
program #20239 as above). The total integration time 
per pixel is estimated to be ~ 2500 sec. 

The fluxes and their uncertainties are estimated from 
PRF photometry on the median-subtracted images as in 
the case of the MIPS70 data. We used a 20" match- 
ing radius. The positional offsets and possibility of ran- 
dom matches are addressed in §2.11. We adopt the lat- 
est value of the MIPS160/im calibration, 41.7MJy/ster 
(Stansberry et al. 2007). The MIPS 160 /xm calibration 
is uncertain to ~ 15%. As mentioned in § 2.5, we adopt a 
color correction of 1.04 for MIPS160/im . An additional 
10% correction is needed, since for small field photom- 
etry filtering results in an underestimated flux densities 
(see MIPS Handbook). We confirm that this is present in 
our data by comparing the pixel values at the positions 
of a few of the brightest sources (including MIPS15977, 
which is detected in both the field map and the pho- 
tometry data). Therefore, all APEX-derived fluxes are 
multiplied by 1.14, which includes color correction, as 
well as the above offset. A total of 9 sources have > 3 u 
detections (1 more has 2.6 cr). Including MIPS22404, this 
adds up to 11 sources (23% of the sample). For the rest, 
we use the 3 a upper limits in all analysis. The typical 
1 (T uncertainty of our MIPS 160 /im data is ^ 9mJy. 

As an illustration of our 160 /xm data. Figure 1 shows 
the MIPS 24 /im images, centered on our target and over- 
laid with the 160 /im contours. We show the 11 sources 
with highest signal-to-noise ratio at 160 /im , where we 
combined the XFLS map (Frayer et al. 2006a) with our 
photometry data. An important point we see in Figure 1 
is that in some cases, the shape of the 160 /im contours 
(e.g. MIPS289, MIPS8184, MIPS22404) suggests that 
more than one source contributes significantly. We ad- 
dress this issue in § 2.10. 

2.7. MAMBO 1.2 mm 

Observations were carried out at the IRAM 30 m tele- 
scope using the 117 element version of the Max Planck 
Millimeter Bolometer (MAMBO) array (Kreysa et al. 
1998) operating at a wavelength of 1.2 mm, during sev- 
eral pool observing periods between October 2004 and 
March 2007. On-off observations were typically obtained 
in blocks of 6 scans of 20 subscans each, and repeated 
in later observing nights until reaching a total rms near 
0.66 mJy or a 5a detection. The median 1.2mm la 
rms noise level reached for the sample is 0.59 mJy (range 
0.35-0.82 mJy). The data were taken during pool pe- 
riods with good atmospheric transmission Ti.2mm < 0.3 
and low or medium skynoise, for a typical total observing 
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Fig. 1. — The 11 160 -detected sources (our IRS sources are 
marked with small red crosses). Here we combine the XFLS map 
with our targeted photometry data. The 24 /xm image is overlaid 
with the MIPS160/im 1 — 5 a" contours. In some cases, multiple 
24 /im sources contribute to the 160 /tm emission (see §2.10 for 
details). The boxes are 120" X 120", with the N-E orientation as 
shown. 



time of 2h per source spread over several periods. Data 
were reduced with standard procedures in the MOPSIC 
package developed by R. Zylka. The flux calibration is 
based on multiple observations of the IRAM flux cali- 
brators during the pool sessions and monitoring of the 
1.2 mm atmospheric transmission through regular sky- 
dips, and is believed to be accurate to ~15%. Correlated 
sky noise was subtracted from the active channel us- 
ing average signals from surrounding channels. All data 
taken during a several weeks period affected by technical 
problems after a change of the telescope control system 
were discarded, as well as individual scans likely sub- 
ject to telescope or instrument problems according to 
the telescope log and inspection of the data themselves. 

MAMBO 1.2 mm observations have been obtained for 
44/48 = 92% of our IRS sample with redshifts. Of these, 
7 are detected at >3ct (14 at >2g). All sources with- 
out individual 1.2 mm detections, have a weighted mean 
signal of (S'l. 2mm) = 0.47 ± 0.09 mJy (the uncertainty is 
error- in-mean, i.e. rms/-\/7V), and are therefore detected 
as a sample at > 5 ct significance. 

Part of these observations were previously presented in 
Lutz et al. (2005). Lutz et al. (2005) obtained MAMBO 
photometry for a sample of bright XFLS 24 /im sources 
with i?(24,0.7) similar to our IRS sample, but without 
the additional i?(24,8) constraint. Including only targets 
from our IRS sample. Tablet lists all available 1.2mm 
fluxes and their statistical uncertainties. 

2.8. YLA 20 cm 

Radio 1.4 GHz (20 cm) fluxes are available from the 
Very Large Array (VLA) map of the XFLS field (Con- 
don et al. 2003). This was done in B-array with a 
FWHM of 4.8" at 1.4 GHz. The typical rms of the 
map is ~23/iJy. We matched our sources to the 3.9cr 
('5'i.4GHz >90/iJy) catalog (J.Condon, private communi- 
cation). We assumed a matching radius of 2.5" (the size 
of the MIPS24/im beam). We find counterparts for 30 
of our sources, where the positional offset with respect to 
the 24 /im positions was in the range 0.2" -1.6" with an 
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average value of 0.9". One additional source, MIPS8342 
is included despite having a positional offset of 3.6" as its 
radio image looks like this might be an extended source 
associated with the 24 /um source. The smaller XFLS 
verification field was also imaged to an rms of ~ 8 ^ Jy us- 
ing the Westerbork Synthesis Radio Telescope (WSRT) 
(Morganti ct al. 2004). Only 3 of our sources (MIPS8245, 
MIPS8268, and MIPS15977) fall in this field. For these, 
we use the deeper observations (MIPS8268 is undetected 
in the shallow map). This brings our radio detections to 
67% of the sample. For the rest, we use 3 a upper limits. 
The 1.4 GHz radio fluxes are listed in Table 1. 

2.9. GMRT 610MHz 

In addition to the 1.4 GHz fluxes, 48% of our sam- 
ple is detected in the Giant Meterwave Radio Tele- 
scope (GMRT) 610 MHz map of the XFLS fleld (Garn 
et al. 2007). The typical 3cr limit of the catalogue is 
~0.3mJy. We describe the radio spectrum as a power 
law of the form F,^ ocz/", where a is the spectral index. 
At the redshifts of our sample, the GMRT data (com- 
bined with the VLA data) allows us to constrain the rest- 
frame 1.4 -5.0 GHz radio spectral indices. The avail- 
able 610 MHz fluxes as well as the radio spectral indices 
are listed in Table 1. Lastly, we would like to note that 
all fluxes (70 /xm -610 MHz) shown in Table 1 represent 
total fluxes, therefore no further aperture matching is 
required. 

2.10. Effect of multiple sources 

Our IRS spectra cover a spatial scale of ~5", com- 
parable to the MIPS24/xm beam. By comparison the 
MAMB01.2mm, MIPS70^m and 160 /xm beams are 
~ 11", 20", and 40" respectively. The large beams, espe- 
cially of MIPS 160 /im , suggest the possibility of multiple 
sources contributing to the measured flux (see also Fig- 
ure 1). In Figure 2, we quantify this by measuring the 
percentage of the total 24 /im flux within a given beam 
that is contributed by the brightest 24 /^m source. This 
analysis has been applied previously in the case of ISO 
170 /xm sources (Sajina et al. 2006), as well as SMGs 
(Pope et al. 2006). 

Given the low signal-to-noise ratio of our sources, 
we assume multiple sources are significant only where 
the percent flux difference is >30%. Within the 
MlPS160/Lim beam, 15 sources meet this criterion (31% 
of the sample). Within the MIPS70/xm beam, 5 sources 
have multiple 24 /im sources (11% of the sample). Fi- 
nally, only 1 of our source has multiples within the 
MAMBO beam (2% of our sample). For sources with- 
out individual detections, attempting to deblend mul- 
tiple sources is futile. For the sources with individual 
160 /im and/or 70 /im detections and multiple 24 /im 
contributors, we discuss our deblending strategy on a 
source- by-source basis in §2.10.1. 

2.10.1. De-blending MIPS70 iim and 160 fim sources 

MIPSllO Here there is an additional 24 /im source 
(0.8 mJy) close to our source (1.8 mJy). However, these 
sources are within 2" of each other (this is the only multi- 
ple source within a MAMBO beam). This close, we can- 
not disentangle the contribution of the two 24 /xm sources 
and hence treat them as one. There is an ofl^set between 
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Fig. 2. — The relevance of multiple 24 /im sources within the 
MIPS70 jum , MIPS160 /^m , and MAMB01.2 mm beams. The x- 

axis is calculated as 100(1 — O'^'' sources usually cor- 

respond to S2A,max- If thc percentage flux dilference is within 
~30%, it is lilcely included in the far-lR flux error, and is negligi- 
ble. Differences ^ 30% are lilcely to have two or more IR sources 
significantly contributing to the measured flux. 



the 70 /im APEX source and our nominal source position 
(7.2"); however, the 70 /im peak does not correspond to 
any other 24 /xm source, and hence is still most likely due 
to our source (see e.g. §2.11). 

MIPS279 - This is a chain of three 24 /xm sources start- 
ing with our target. They arc separated by ^ 10" be- 
tween the first and second and second and third sources. 
The 70 /im source is centered on our target (positional 
offset 1.0"), and hence we assign the full 70 /xm flux to 
our source. Of the three, only our source has a 20 cm ra- 
dio detection. The radio flux is 6 cr the local map rms sug- 
gesting that the undetected radio emission of the other 
sources is significantly lower than our target's. Therefore 
most likely the MIPS160/xm far-IR emission is domi- 
nated by our target as well. 

MIPS289 Our source appears blended with a nearby 
bright 54 ± lOmJy source (~ 20" away). We attempt to 
deblend the two by combining two 2D Gaussians, and 
find that our source has about half the flux of its neigh- 
bor. We therefore assign it a flux of ~ 27 ± 10 mJy. 
MIPS8184 - This source is part of a group of 6 
MIPS24/im sources ranging from 5*24 = 0.4 1.2 mJy 
within the 160 /im beam. Our source position however 
is the strongest 70 /xm source. The ratio between the 
peak 70 /xm flux at the position of our source and the 
next brightest peak is 1.5. Therefore we assume that our 
source dominates the 160 /xm emission as well. 
MIPS8196 - Here there are three 24 /im sources, one 
5.2" from our target and the other a further 10" from 
the second source. The 70 /im emission is centered in 
between and seems to be largely contributed to by the 
later two sources. Therefore we treat this source as a 
70 /im non-detection. 

MIPS8342 - This system has a nearby (~ 20") neighbor 
of roughly equal 24 /im flux. The centroid of the 160 /im 

emission lies roughly in-between the two 24 /im sources 
(see Figure 1). The ratio of the 70 /im flux ratio between 
the position of our source and the secondary source is 1.2, 
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which is also consistent with the two sources having com- 
parable far-IR emission. It seems reasonable to assume 
that our source contributes half of the total 160 /um flux 
(29±llmJy). We therefore take its deblended 160 /xm 
flux to be half of that, or 14~ 11 mJy. 
MIPS22404 At 160 /xm this is a large blob encompass- 
ing at least four 24 /im sources (sec Figure 1). Two of 
the sources are far enough away that they are excluded 
from the PRF-derived 160 /zm flux. However, the near- 
est source to our target is 7.2" away, with the 70 /xm 
emission centered roughly in between our source and its 
neighbour. Moreover, both are radio sources. It is likely 
that both sources contribute to the far-IR emission, al- 
though their relative contribution is uncertain. For sim- 
plicity, we split both the 70i^m and 160 /xm fluxes be- 
tween our target and its neighbor. 

Due to the uncertain nature of the above deblending 
procedure, we prefer to quote the original APEX fluxes 
in Table 1, but use the deblended values in all subsequent 
analysis. 

2.11. Random matches 

In § 2.5 and § 2.6, we discussed the flux determination 
for the MIPS70/im and MIPS160^m observations re- 
spectively. In both cases, we ran APEX first to find all 
sources in the images and then matched those against 
our 24 /im source positions. Table 1 shows the fluxes of 
the detected sources as well as the positional offsets be- 
tween the APEX-derived 70 and 160 /im sources and the 
reference 24 /xm positions. In this section, we address 
the chances of random associations in this procedure for 
cases where there is a single 24 /im source. The cases, 
where multiple 24 /im sources may contribute to the ob- 
served far-IR emission are discussed in §2.10.1. 

We apply the P-statistic (Downes et al. 1986), which 
estimates the random chance of one or more sources be- 
ing found within a radius of a given source. It is de- 
flned as P = 1 - exp^—irnO^), where n is the source den- 
sity above a given flux level, and 6 is the search ra- 
dius. Neglecting the cases where multiple 24 fxm sources 
can contribute to the 70 /zm flux we find that the posi- 
tional offsets vary between 0.1" and 6.9", with an aver- 
age value of 2.4". This is comparable to the nominal 
MIPS70/xm positional uncertainty which is 2.6" (1 cr; 
Prayer et al. 2006a). We compute the P-statistic for 
the source, which has the highest chance to be a ran- 
dom match. This is MIPS22600 where a 6.4 mJy 70 /im 
source is found at 6.9" from the 24 /im source. We use 
the 70 /im counts from Frayer et al. (2006b) which give 
2.6 X lO^gals/ster > 6.68 mJy (with a counts bin ranging 
from 6.0 to 7.5 mJy). Given 6' = 6.9" for MIPS22600, we 
find P = 0.009. Given that the faintest detected source 
(MIPS8242 with 4.8 mJy) is already 75% the flux of 
MIPS22600, but with smaller positional offset, the rest 
of the sample would have even lower probabilities for be- 
ing random associations. The low P values mean that 
our 70 /im detections are not random. 

For MIPS160/im we have 11 detections (sec §2.6); 
however. 5 of these arc blended sources (sec §2.10.1) 
making this statistic difficult to apply. For the remaining 
6 sources, we find that the positional offset between the 
160 /im and 24 /im source ranges from 1.5" (MIPS22482) 
to 9.8" (MIPS22651), with an average offset of 6.3". This 
is comparable to the typical positional accuracy of 5.2" 



at 160/im(lcr; Frayer et al. 2006a). MIPS15880 repre- 
sents the 'least likely' detection having among the great- 
est separation between the 24 /zm and 160 /im sources 
(9.3") as weU as the lowest flux (23.8 mJy). For 160 /xm, 
there are no published source counts that go down to our 
lowest detected fluxes (~30mJy). Based on the model 
counts from Dole et al. (2004a), the range of acceptable 
models gives ~ 1 3 x lO^gals/ster having ^ 30mJy. For 
MIPS15880, a 23.8 mJy source within 9.3", we find that 
the above range gives P = 0.006-0.019. The rest of the 
non-blended 160 /zm detections, would have P -values no 
worse than the above. Therefore, we conclude that likely 
none of the 160 /im sources are random associations of 
our 24 /im sources. 

2.12. Stacking 70fi.m undetected sources 

Of the weak-PAH, z> 1.5 sources nearly half (12/27) 
are undetected at 70 /im (see Table 2). The aver- 
age of the aperture photometry measurements for these 
sources (see Paper II) is ^ 2 mJy which is about four 
times lower than the average of the 70 /im -detected weak- 
PAH sources at 2;>1.5 (-^lOmJy). For another mea- 
surement of the average flux of these sources we also 
stack the 70 /zm image at the positions of these sources. 
Our stacking largely follows the technique of Huynh et 
al. (2007b). Here we summarize the key points. Be- 
fore stacking, we subtract all bright (> 5 a) sources from 
the median-subtracted image. To estimate the reliabil- 
ity of our stacking result, we also perform 500 stacks of 
24 random position. We find a stacked flux of 1.95 mJy. 
However, the standard deviation of the offset stacks is 
0.68 mJy making this only a 2.8 cr detection. 

2.13. Stacking 160 jim undetected sources 

Table 2 shows that only 3 of the 27 z>1.5 som'ces 
are individually detected at 160 /zm. These are all also 
detected at 70 /zm. There are a total of 24 ^;>1.5, 
weak-PAH sources without individual 160 /zm detections. 
However, we note that our 160 /im -detected sources have 
fluxes already roughly consistent with the model con- 
fusion limit of 160 /zm (Dole et al. 2004). This sug- 
gests that stacking is unlikely to improve this signifi- 
cantly. Nevertheless, using the same stacking technique 
as describe in §2.12, we arrive at an average 160 /zm 
flux of 9.2 mJy. The standard deviation of the offset 
stacks is 5.9 mJy. Therefore these sources are unde- 
tected in the stacked image. We adopt the 2 a upper 
limit for these sources (i.e. 12mJy). It is likely that 
the 70 /zm -undetected among these sources have lower 
160 /im -fluxes than the 70 /zm -detected ones; however, 
we flnd that we cannot distinguish between these sub- 
samples in our stacking for the reasons described above. 

3. RESULTS 

3.1. Overview of our AGN/starburst diagnostics 

There are two distinct questions we wish to address 
here: 1) whether or not our sources contain AGN and/or 
starbursts, and 2) what is the relative contribution of 
the two processes to their bolometric power output. The 
presence of AGN in the bulk of our sample was origi- 
nally asserted on the basis of their typically low PAH 
equivalent widths, and high 6 /tm hot dust continuum 
(see Paperll). Conversely, the strong PAH emission in a 
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quarter of the sample, and the presence of PAH in the 
stacked spectra of the weak-PAH sources suggest that 
star-formation is hkely present in the majority of our 
sources (Paperll). The presence of AGN in many of the 
weak-PAH sources was independently shown on the ba- 
sis of their AGN-like radio luminosities and even, in two 
cases, radio jets (Sajina et al. 2007b). 

In this section we further test for the presence of star- 
bursts by the detectability of cold (< 50 K) dust emission 
typically associated with this process. We further ex- 
pand the discussion on the presence of radio AGN on 
the basis of radio luminosities, spectral indices, and the 
ir-radio correlation. For a sub-set of the sample with 
available optical/near-IR spectra, we present optical line 
diagnostics which could reveal the presence of AGN. 

Finally, we use SED decomposition to estimate the rel- 
ative contribution of AGN and starbursts to Lir . To do 
so, we make the simplifying assumption that the cold 
dust continuum is entirely due to star-formation while 
the warm/hot dust continuum is entirely due to AGN. 
The derived starburst component luminosities are the 
sum of the cold dust and PAH luminosities. We com- 
pare our results from the above analysis with the SEDs 
of well known sources and models in order to assess the 
level of uncertainty in the relative AGN/starburst frac- 
tions of LiR. In particular, we address the questions: how 
much starburst activity is there in the weak-PAH sources 
and how much AGN activity is there is the strong-PAH 
sources? 

3.2. Far-IR continuum: starburst or obscuration? 

Typically, the far-IR emission is assumed to arise from 
star-formation even in quasars (e.g. Beelen et al. 2006; 
Marshall et al. 2007; PoUetta et al. 2008). However, for 
highly obscured AGN as is the case with our sources (see 
Paperll), there exists the possibility of dust reprocessing 
the typically hot/ warm SED of an AGN into a much 
cooler starburst-like SED (e.g. Levenson et al. 2007). 
Therefore, before we proceed, we look for trends between 
the far-IR continuum emission and the PAH equivalent 
widths (indicative of starburst activity) and 9.7 /im sili- 
cate feature optical depth (rg.y) indicative of obscuration. 

As shown in Table 2, both the 160 /zm and 1.2 mm de- 
tectability of our sample are rather low (21% and 16% 
respectively). However, in both cases, the detectability is 
higher for the strong-PAH vs. weak-PAH z > 1.5 sources 
(e.g. 50% vs. 11% for the 1.2mm data). This is shown 
in Figure 3left. To assess the strength of the perceived 
trend, we perform a Spearman rank correlation. We find 
p =0.569 with probability of correlation of 99.9%. By 
contrast, the MAMBO 1.2 mm flux vs. the optical depth 
of the 9.7 fim silicate feature (Figure 3right) has a Spear- 
man rank correlation of p — 0.146 with probability of cor- 
relation of 68.2%. This suggests that there is a significant 
trend of strong-PAH sources to be strong cold-dust emit- 
ters, but no such trends exists with respect to the degree 
of obscuration as quantified by the depth of the silicate 
feature. This conclusion supports the finding of Desai et 
al. (2007) who show that in local ULIRGs, the degree of 
obscuration and the far-IR colors are not related. 

3.3. Presence of AGN: Lradio md spectral indices 

Table 3 gives the rest- frame 1.4 GHz luminosities of 
our sample, where the k - correction is based on their 
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Fig. 3.— Left: MAMBO 1.2 mm flux vs. EW7.7. To minimize 
the effects of negative fe-correction, we only plot 2>1.5 sources. 
Right: MAMBO 1.2 mm flux vs. T9.7. We have added 0.1 to all r 
values in order to be able to plot the sources with rg.7 = 0.0. In all 
cases, the red symbols are strong-PAH sources {EW7.7 >0.8/^m), 
while the blue symbols are weak-PAH sources (EW7.7 < 0.8 /im). 



observed a'j'^^Q^^ spectral indices (see Table 1), or 
assuming ai"'^4GH^^ =-0.7 for the sources not detected 
in 610 MHz.' In Sajina et al. (2007b), we adopted 
Liaghz > lO^^W/Hz as our radio-loud criterion. Fig- 
ure 4 shows the radio luminosities for the entire sam- 
ple. It is obvious that most of the non-radio loud, radio 
detected sources still have significant radio luminosities 
(ii.4GHz ^ lO^'' - 10^5 W/Hz). Most sources in this range 
are likely AGN powered as well. For example, the radio 
luminosities of SMGs, which were found to be extended 
(i.e. starburst-powered) in high-resolution radio imaging 
all have radio luminosities < 10^"* W/Hz (Chapman et al. 
2004). 
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Fig. 4. — Distribution of radio luminosities vs. redshift for 
our sample. The dotted line shows the radio-loud limit we adopt 
(Sajina et al. 2007b). The dashed line shows the XFLS 4a limit of 
90 /ijy (assuming Oi^^^cfuz =-0.7). The upper limits are our radio- 
undetected sources which by doflnition lie on the dashed line. 

The radio spectral index, a, is typically assumed to 
be ~-0.7 for star-forming galaxies (Condon 1992), and 
^-0.5 for quasars (Stern et al. 2000). Recent blank-sky 
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radio observations are consistent this range (Bondi et 
al. 2007; Huynh et al. 2007a). Steep spectra (a~-l) 
are however the norm for high-z radio galaxies (HzRGs). 
Flat or inverted spectra (essentially a > -0.4) are strongly 
indicative of AGN. 

Figure 5 shows the distribution of af^" for all sources 
with both VLA and GMRT detections, as well as the 
radio-loud and non-radio-loud sources separately. The 
bulk of our radio-detected, but not radio-loud sources 
show fairly standard spectral indices (a'j'-'^ ^--^-0.8). As 
discussed in Sajina et al. (2007b), our radio-loud sample 
tend to have steep radio spectral indices {af^^ ^-1). This 
is consistent with the steep spectra found for a sample of 
z ^ 2 Type2 quasar candidates (Martfnez-Sansigre et al. 
2006). Magliocchetti et al. (2008) also find that the bulk 
of their optically- faint 24 /xm -bright sources show steep 
spectra (a^^^^ ^ -1.0)). 
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Fig. 5. — Distribution of radio spectral index for the sources with 
both VLA 1.4 GHz and GMRT 610 MHz detections. We also show 
the radio-loud and non-radio-loud sources separately. 

Three of our sources (MIPS227, MIPS16059i°, and 
MIPS22303) are consistent with flat spectra (a^i4°|g^^ > - 
0.4), by virtue of being undetected at 610 MHz. All of 
these sources are weak-PAH sources. For convenience, 
we have added a column in Table 3 indicating whether or 
not a source shows signs of the likely presence of a radio 
AGN (i.e. is it radio- loud and/or has a flat spectrum). 
Given that the typical uncertainty in the radio luminos- 
ity is ^0.2dex, we also identified borderline radio-loud 
sources (i.e. those with log Li.4ghz = 24.8-25.0). As dis- 
cussed above, these luminosities are also more than a 
factor of five higher than the strongest known starburst- 
powered radio luminosities. These are also very likely 
AGN-powered. 

We find that 12/34 (35%) of our weak-PAH sources 
show signs of AGN. Of the strong-PAH sources, 4/14 

This source shows an inverted spectrum, given its 1.4 GHz flux 
of 0.57 mjy and no 610 MHz detection. The formal 610 MHz detec- 
tion limit (3 (T we adopt is 0.3 mJy (see § 2.9). However, MIPS16059 
is in a noisier part of the GMRT map and hence we estimate that 
a detection limit of ~0.6mjy is more appropriate, which moans 
the radio spectrum of MIPS16059 may also be flat. 



(29%) show signs of radio AGN. This shows that there 
is no clear trend between the PAH strength and radio- 
AGN detectability. However, in Sajina et al. (2007b), 
we showed that of the high-z weak-PAH sources, only 
sources with T9.7 > 1 showed radio-loud luminosities 
(specifically 40% of the sources in this category were 
radio- loud). The reason for this finding is still unclear. 

3.4. Presence of AGN: optical line diagnostics 
3.4.1. High ionization lines 

Figure 6 shows the 13 available optical (rest-frame UV) 
spectra for our sample. These include 12 Keck LRIS and 
DEIMOS spectra (§ 2.2). . The rest-frame UV spectrum 
of MIPS22204^^ was presented in Martmez-Sansigre et 
al. (2008). For completeness, we include it here as well. 

MIPSllO shows the [Nev] 3426A fine, which re- 
quires an AGN for excitation. MIPS227, MIPS15928, 
MIPS16080, and MIPS22204 aU show the Civl550A 
emission line. Our velocity resolution at 1550Ais 
'--^280km/s. In all cases, the line is resolved with 
a FWHM of -1241km/s, 2058km/s, 2507km/s, 
and 1305 km/s respectively for MIPS227, MIPS15928, 
MIPS16080, and MIPS22204. These widths are 
consistent with the Narrow-Line Region (NLR) of 
AGN(Sulentic & Marziani 1999). All these sources 
are therefore consistent with hosting obscured AGN. 
With the exception of MIPS15928 (which is a border- 
line strong-PAH source, see Paperll), all these sources 
are have weak-PAH mid-IR spectra. By contrast, the 
five sources which do not show signs of AGN (MIPS283, 
MIPS8207, MIPS22404, MIPS22554, and MIPS22600) 
are all low-z, strong-PAH sources. The spectrum of 
MIPS133 is too red for the above diagnostic lines. How- 
ever, we observe the [Oiii] line (see §3.4.2). MIPS22699 
shows no lines in its DEIMOS spectrum, which is not sur- 
prising if we assume its IRS-derived redshift of z = 2.6. 
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This source is the same as AMS13 in Martfnez-Sansigre et al. 
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Fig. 6. — The available rest-frame UV/optical spectra for our 
sample. Where detected, a continuum has been subtracted. The 
y-axis are in arbitrary units. The spectra for MIPSllO, MIPS133, 
and MIPS15928 have been smoothed by a factor of 5 for clarity. 



Figure 6 — Continued. 

3.4.2. Optical line ratios 

Figure 7 shows the rest-frame optical spectra for the 
6 sources^^ with Keck NIRSPEC data and Figure 8 
shows the same for 5 sources with Gemini NIRI data. 
The H/3, [Oiii]AA4959,5007A as well as the Ha6365A 
and [Nii]AA6549,6583A complex are all marked. These 
emission lines are fit with Gaussian profiles of fixed cen- 
tral wavelength but variable width and amplitude. The 
complex consisting of Ha and the two [Nil] lines is fit 
simultaneously with three Gaussians all of fixed central 
wavelength, and includes a single linear continuum com- 
ponent. We determine the redshifts to fit the best both 
the H - band and K - band lines and estimate the near- 
IR redshift uncertainties at ±0.001. The only exception 
is MIPS15949 where we found different redshifts were 
needed for the two (see §3.4.3). In this fitting we use 
the IDL routine MPFIT. The principle line strengths and 
line ratios are given in Table 4. The FWHM values have 
been corrected for the respective instrument resolutions 
(see §2.3). 

We adopt the following criteria for detecting an optical 
AGN: 

1) FWHM(Hq;) >2000km/s (e.g. Zakamska et al. 2003, 
and references therein) 

2) log([Nii]/Ha) > -0.2. - This is based on Kauffmann et 
al. (2003), which cut likely includes composite starburst- 
AGN sources (see e.g. Kewley et al. 2006). This criterion 

12 An additional source observed with NIRSPEC was MIPS78; 
however, no lines were detected in its spectrum. Its IRS-based 
redshift is 2: =2.65, which means that Ha falls outside its observed 
range. It is also quite faint {nifj =21.98; Dasyra et al. 2008). 



is used only when log([Oiii]/H/3) is not available. 
3) When both ratios are available, we use the Kewley et 
al. (2006) criteria to determine whether or not a source 
is a starburst, an AGN, or a composite source. 

Table 4 shows the classification of our sources accord- 
ing to the above criteria. We find that only 1 of our 
11 sources (MIPS15949) can be classified as a broad-line 
(FWHM> 20G0km/s) AGN. Another potentially broad- 
line source is MIPS22204; however, we suspect it to be a 
blend of two sources nearly overlapping on the sky and 
with a small velocity separation (Aw ~ 500km/s). The 
redshift (z — 1.974) is based on the Lya and CiV detec- 
tion presented in Martfnez-Sansigre et al. (2008). Based 
on its Civ line we classified it as an AGN in § 3.4.1. 
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Fig. 7.— The Keck NIRSPEC data. The y-axis units are 
lO""*^* erg/s/cm^/A. We simultaneously fit the principle lines with 
Gaussian profiles of fixed central wavelength (see text for details). 
The best-fit linear continua are shown by the blue lines, while the 
red lines show the best-fit continuum-|-emission lines model. The 
fit for MIPS22204 is omitted as it appears to be a blended source. 
For each source, the top panel shows the 1-D spectra, while the 
bottom shows the 2-D spectrum. The later are difference images 
(to subtract the sky lines) and are based on individual observa- 
tions, not their mean as the 1-D spectra. The 2-D spectra have 
been smoothed with a 5 pixel kernel to highlight the features. 

We find that, for the sources with optical/near-IR spec- 
tra, 3/9 (33%) of the strong- PAH sources show signs of an 
optical AGN, while 9/13 (69%) of the weak-PAH sources 
show signs of AGN. The much higher fraction of the later 
indicates that there is at least some correspondence be- 
tween the mid-IR AGN diagnostics (i.e. PAH strength) 
and the optical AGN diagnostics. We cannot fully ex- 
clude the possibility that some of even our weak-PAH 
sources are pure starbursts with unusually strong mid-IR 
continua (e.g. MIPS429, MIPS8196); however, the bulk 
of the weak-PAH sources do show evidence of obscured 
AGN. 

3.4.3. Asymmetric [Olll] profiles 

An additional, diagnostic is the profile shape of the 
[OIII] 5OO7A line. Because of the poor signal to noise 
in many of the spectra, this is difficult to ascertain in 
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Fig. 9.— The [OIII] profile of MIPS16080 showing a broader 
blueshifted component in addition to the narrow lines. 



Fig. 8. — The same as Figure 7, but for the Gemini NIRI spectra. 
The only exception is MIPS8327, whose iC-band spectrum was 
taken with NIRSPEC. 



most sources; however, some sources show clearly the 
blue-asymmetry frequently observed in the [OIII] lines 
of AGN (Heckman et al. 1981; Rice et al. 2006). This 
asymmetric profile is believed to be a combination of 
excess blue wings due to outflows, and obscuration of 
the red wings (Heckman et al. 1981). Figure 9 shows 
a close-up of the [Oiii] line of MIPS16080, which shows 
this asymmetry most clearly. We have performed a si- 
multaneous fit of 3 Gaussians two of which are fixed to 
the [Oiii] 5G07A and 4959A lines, the central wavelength 
of the third line is not fixed but is restricted to be in 
between the other two. The best-fit additional broad 
component has FWHM = 1849km/s and a velocity off- 
set from the 5007A line of Av = -261 km/s. 

Another potential case of such an outflow is 
MIPS15949. We find that there the [Olll] redshift is 
2.116, while the Ha redshift is 2.122. This blueshift of 
the [Oiii] line, if attributed to an outflow, imphes outflow 
speed of 576 km/s. 

3.4.4. Emission-line based AGN luminosities 

Two effects make the determination of the intrinsic 
AGN L[oiii] uncertain in our typically composite sources: 
the unknown level of star-formation contribution to this 
line, and the unknown level of extinction. Despite these 
caveats, the observed [Oiii] luminosities of our sources 
can still serve as a rough gauge of the strength of the 
central source. The [OiiiJ luminosities are in the range 
^;^0*-56-9.02l^^ or - 10''2-2-42.6g^.g/g ^i^ese are com- 
parable to the more luminous of the Type 2 quasar can- 
didates selected in the SDSS by Zakamska et al. (2003). 
Without reddening correction on the emission lines, Net- 
zer et al. (2006) derive a relationship between the [Oiii] 
and X-ray luminosities of Type2 quasars: log(L[oiii] 



/ L2-iofeey) = (16.5±2.9) - (0.42±0.07) * L2-iokeV- 
Based on this, we expect intrinsic X-ray luminosities in 
the range ~ 10'*'''^~^^ °erg/s for our sources with [Oiii] 
detections. Even accounting for the flux calibration un- 
certainty in the spectra (see §2.4), our sources are con- 
sistent with quasar-strength AGN. The Chandra X-ray 
observations of our sample and obscuring columns are 
presented in Bauer et al. (2008, in prep.). 

3.4.5. Optical Extinction and SFR 

We have four sources with an estimate of the Ha /HP 
ratio. These are: MIPS8196, MIPS8327, MIPS15949, 
and MIPS16059 with ratios of 2.5, 2.0, 9.8, and > 25.7 
respectively. In a comprehensive study of 2000 Syl AGN, 
Zhou et al. (2006) find a mean value of Ha/H/3 = 3.03, 
consistent with theoretical expectations, (see Osterbrock 
1989). The observed range in Zhou et al. (2006) is 
roughly Ha/H/3~2-5. Our measured ratios may suf- 
fer some systematic uncertainty due to the uncertainties 
in the absolute fiux calibration of the H and K bands 
(see § 2.4). We can however conclude that MIPS8196 and 
MIPS8327 have low levels of optical extinction (if any), 
while MIPS15949 and MIPS16059 show significant ex- 
tinction. Using a standard Milky Way extinction curve, 
we obtain E(B-V) = 1.2 and E(B-V) > 2.2 for MIPS15949 
and MIPS16059 respectively. For comparison. Brand et 
al. (2007) find that all their IR-bright, optically-faint 
z ^^2 sources with H/3 coverage, show strong extinction 
in the range E(B-V) = 1.0-1.9. 

The two sources without optical extinction, MIPS8196 
and MIPS8327, have respectively rg.v = 1.3 and 2.4 (Pa- 
perll). The two sources with strong optical extinction, 
MIPS15949 and MIPS16059, have respectively T9.7 = 0.0 
and 2.7. This suggest that there is no correspondence be- 
tween the optical extinction values derived above and the 
strength of the 9.7 /im silicate absorption feature. This 
is consistent with the conclusions of Brand et al. (2007). 
In both cases, we are dealing with fairly small number 
statistics; however, this conclusion is not surprising given 
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the expected comphcated dust/power source geometries 
in these sources. 

The average Lna of our sample (excluding the broad- 
line MIPS15949) is 2x IG'^'^eTg/s (Table 4). This com- 
pares well with the non-broad line sources in Brand 
et al. (2007). Based on Kennicutt (1998), we derive 
(SFR(Ha)) =9.8Mo/yr (range 1.3-42.4M0/yr). This 
is a lower limit given the expected extinction in most 
of these sources. For MIPS8196 and MIPS8327 (two 
weak-PAH sources without significant optical extinc- 
tion), we find SFRs of 22MQ/yr and 26MQ/yr respec- 
tively. These modest SFRs are consistent with the far-IR 
non-detection of both these sources (see Table 1). We re- 
turn to the question of IR-derived SFRs in the following 
section. 

3.5. Far-IR SED fitting 

In §3.2, we show that, especially in strong-PAH 
sources, but also in some weak-PAH sources, substan- 
tial cold dust emission is present further supporting the 
presence of starbursts in these sources. Here, we fit the 
available far-IR data (70 /im, 160 /im, and 1.2 mm) in 
order to both determine the total infrared luminosities 
of our sources, and constrain the relative contribution of 
AGN and starburst processes to those luminosities. The 
far-IR SED of infrared-luminous galaxies is frequently 
modeled as a composite of two or more different tem- 
perature modified blackbody components (e.g. Dunne & 
Eales 2001; Klaas et al. 2001). The modification is given 
by . The emissivity index, /3, ranges between ~1- 
2. The greybody approach is only an approximation to 
a much more complex radiative transfer problem, but 
is useful in characterizing far-IR SEDs when only a few 
data points are available. We only have data at 70 /im , 
160 /im, and 1200 ^m, which, at 2^2, correspond to 
Arost ~ 20, 50 and 400 ^m . 

In Paperll, we already discussed that our mid-IR con- 
tinuum can be thought of as a sum of different tempera- 
ture dust components ranging from ~ 1500 K to ~ 100 K. 
Here, we assume that this mid-IR continuum is entirely 
AGN-powered, and use this model and set the cold- 
est AGN component at Touter = 100K, /3=1.5. This 
choice gives good fits for the vast majority (83%) of 
the sources; however, some exceptions exist. In par- 
ticular, for MIPS506, MIPS8242, and MIPS15840 we 
used Toutcr = 150K, while for MIPS8342, MIPS15928, 
MIPS16080, MIPS16095, and MIPS22277 we found that 
Toutcr = 80 K gave better agreement with the data. This 
continuum of temperatures was additionally subjected to 
the extinction of a level determined by the 9.7 /xm sili- 
cate feature depth (see Paperll). Our fitting procedure 
in Paperll includes the uncertainty in the determination 
of this mid-IR continuum which was especially strong for 
high-z, strong-PAH sources. We find however, that the 
uncertainty in the total Lagn requires a further consider- 
ation of the uncertainty in Toutor, which is the dominant 
component for steep-slope sources. This is difficult to 
quantify robustly, but we find that a range of ± 10 K in 
its value roughly corresponds to the uncertainty in the 
MIPS70/im points (rest-frame ~20/xm). In Figure 10, 
we show the AGN component with its 1 a uncertainty, 
which includes both the fitting uncertainty from Paperll 
and this reasonable spread in the value of Toutor- We 
further discuss the AGN-powered mid-IR SEDs in §3.7. 



Before fitting for the starburst-powered cold dust emis- 
sion, we subtract the AGN-component as described 
above from the MIPS70 and MIPS160/im points. Be- 
cause of the typical far-IR non-detection of our sources, 
or detection of 1 - 2 points only, we cannot justify fitting 
more than a single greybody component (there are three 
exceptions here as described in Figure 10). We use the 
IDL routine MPFIT to fit for this cold dust component and 
determine its 1 a spread. In the fit, the amplitude, and 
temperature are free parameters while we fix the emis- 
sivity index, (3 = 2. Degeneracy between and /? (see 
e.g. Sajina et al. 2006) means that this high /3 value leads 
to cooler temperatures. 




Fig. 10.— The ~ 1-1000 /xm SEDs of our sources with far-IR 
detections. The blue curve shows the best-fit to the mid-IR spec- 
tra (see Paperll). The dotted red curve shows the best- fit cold dust 
component with the pink shaded region being the 1 a scatter. In a 
few cases with all three far-IR points detected, a second greybody 
component was included to give a 'smoother' SED, but without sig- 
nificantly contributing to Lir. For MIPS289 and MIPS22651 this 
component has =60K, while for MIPS15977 it has = 20K. 
Upper limits are 2 a. In all cases, the solid red curve is the total 
SED. 

Figure 10, shows the best-fit SEDs for all sources where 
at least 2 of the 3 photometric points are detections. For 
the rest, the far-IR is much less constrained. However, 
in the vast majority of these cases, the AGN-component 
dominates (and hence Lir is determined reasonably well 
despite the far-IR non-detections). The average SED of 
the weak-PAH, far-IR non-detections is discussed in § 3.7. 
The integrated 3- 1000 /im luminosities (Lir), as well as 
the luminosities of the AGN and starburst components 
separately are presented in Table 3. The starburst lu- 
minosities listed in Table 3 are the sum of the cold dust 
emission and the PAH emission, determined from our 
mid-IR fits in Paperll. 

3.5.1. Star-formation rates 

Lastly, using the standard relation from Kennicutt 
(1998), we derive the IR-based SFRs for the 160 /im- 
detected sources and upper limits for the majority of the 
sources which are non-detections. The z< 1.5, 160 /im- 
detected sources are typically strong-PAH sources whose 
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implied SFR is - 100 SOOMe/yr. The z > 1.5 160 Aim- 
detected sources typically have SFRs of ^QOOMo/yr. 
We believe that the upper limits for our strong-PAH 
sources are likely close to their SFRs. For the weak-PAH 
sources, while the upper limits allow for significant SFRs 
(up to ~800Mo/yr), in at least two cases we find more 
modest SFRs of ^ 20 - 30 M0 / yr from their Ha luminosi- 
ties and low optical extinction levels (see §3.4.5). This 
is likely to be roughly the lower limit on the SFRs of the 
other weak-PAH non 160 /xm -detected z> 1.5 sources. 

MIPS22530 is the only one source with both 160 /xm 
detection and an Ha luminosity. For this source, 
the Ha derived SFR is 7.3M0/yr, while its IR-bascd 
SFR is 990±400MQ/yr. This implies an extinction of 
A(Hq) = 5.3 (or 4.8 if we take the lower limit on SFRir). 
However, as can be seen in Figure 7, the removal of sky 
lines in this source is somewhat suspect and we might be 
missing some of the Ha flux. Nevertheless, reconciling 
the two SFR estimates clearly requires extreme levels of 
obscuration. For comparison, Brand et al. (2007) find 
A(Ha) > 2.4, 3.8, and 4.6 in three of their ~ 2 ULIRGs 
with limits on the Ha/H/3 ratio. Therefore, although ex- 
treme, the required Ha extinction is not unprecedented 
in similar populations. 

3.6. Lmir -io-LiR 

An important parameter derived from the spectral fit- 
ting is the integrated 3-1000/xm infrared hmiinosity, 
LiR (see Table 3). We find that our z~ 1 sources have 
LiR ^ 10^^ L0, while the bulk of the sample which lies at 
2~2, have Lir^7x 10^^ Lq. The greater luminosities 
of our high-z soiirces arc the usual consequence of a flux- 
limited sample (see Paper I for details of the Malmquist 
bias of our sample) . But how reliable are our Ljr values 
given that we fit the far-IR SED from only three points: 
70 iim , 160 /im , and 1.2 mm? At 2; ~ 1, the 160 /im point 
samples near the peak of the SED distribution, while for 
weak-PAH sources without far-IR detections, the mid-IR 
continuum dominates the infrared luminosity in any case. 
For z^2 MIPS160/^m -detected sources we are missing 
the important ~ 100 fim regime. In practice, our fitting 
procedure typically leads to rather cold dust tempera- 
tures which peak at ~ 100 /im (see Tabic 3 and Fig- 
ure 10). We estimate that if we assume warmer SED, 
i.e. ones that peak at the rest-frame of the MIPS160/Ltm 
observations, then the quoted Lir will be overestimated 
by ^ 0. 1 dex. This is comparable to our ijR uncertainties 
(see Table 3). 

It is well known that the mid-IR luminosity of star- 
forming galaxies is correlated with the total infrared lu- 
minosity, although the exact relation derived is some- 
what affected by selection biases (Chary & Elbaz 2001; 
Takeuchi et al. 2005; Bavouzet et al. 2008). However, 
these relations are optimized for the typical star-forming 
galaxy (Ljr ~ 10^° - 10^^ Lq). It is not obvious how well 
they perform at predicting the infrared luminosities of 
ULIRGs, especially those with significant AGN contri- 
bution. The availability of IRS spectra and far-IR-to- 
mm data allows us to test the conversion between mid- 
IR broadband observations and Ljr for such extreme 
sources. 

In Figure 11, we present the conversion between the 
monochromatic 14 /xm luminosity and Lir as well as 
the rest-frame IRAC8 and MIPS24 luminosities for 



comparison with broadband data. Below we present 
the best-fit linear fits in the log- log plane for the 8 fim , 
14 /xm , and 24 /um luminosities respectively: 

log(LiR) = (2.79 ± 0.36) + (0.83 ± 0.03) log(z/L,,8^m()l) 

log(Lm) = (4.43 ± 0.43) + (0.68 ± 0.04) log(t.i,,i4^„Ce) 
log(LiR) = (4.13 ± 0.54) + (0.71 ± 0.04) log(i.L,,24^m03) 

where the rms spread is respectively 0.11, 0.15, and 
0.17 dex for the three relations. We stress that the above 
are fits to our (extremely biased) sample. They are not 
meant to be all-purpose conversion relations, and should 
only be used for sources of comparable mid-IR lumi- 
nosities to our sources. Another major caveat is that 
for many of these sources, the infrared luminosities are 
upper limits only. Figure 11 shows the expected range 
in Lir for these sources. Accounting for the range in 
Lir, we find that the our sources roughly agree with the 
Bavouzet et al. (2008) 8 /xm conversion relation. How- 
ever, our strong-PAH sources, with 160 /um detections 
all have systematically higher Lir per 8 /xm luminosity. 
For some local ULIRGs the gap is even more extreme. 
This result is not surprising given that lower luminosity 
starbursts have higher PAH contribution to Lir (from a 
few percent up to ~ 20% Smith et al. (2007)), while PAH 
contribute from ^1% up to a few percent of the total IR 
luminosity of ULIRGs (Armus et al. 2007). 

For Li4 vs. Lir (Figure llmiddle), we find the slope is 
significantly different from unity, implying a higher MIR- 
to-FIR ratio at the highest Imninosities. Local conver- 
sion relations for star-forming galaxies suggest a slope ~ 1 
(Chary & Elbaz 2001). As expected, our much shallower 
slope is dominated by the weak-PAH sources with their 
strong mid-IR continua. The strong-PAH sources, for 
the most part, agree reasonably well with the Chary & 
Elbaz (2001) relation. 

Finally, Figure llbottom shows the rest-frame 
MIPS24 /im , whose best-fit relation gives the best agree- 
ment with the lower luminosity conversion relation. The 
same caveats as discussed above for the IRAC8 /xm rela- 
tion apply here as well. 

3.7. The starburst contribution in weak-PAH sources 

Apart from obtaining the total infrared luminosities 
of our sources, the SED fitting described in §3.5 pro- 
vides us with the relative contributions of AGN and star- 
formation to Lir. We wish to compare our empirical 
modeling with other approaches in order to assess the 
level of uncertainty in the Lagn and/or Lsb values we 
obtain. In this section, we address this question for the 
bulk of our sample which consists of z > 1.5, weak-PAH 
sources. We discuss the validity of our decomposition for 
strong-PAH sources in the following section (§ 3.8). 

Figure 12, shows the average SED of the 160 /xm- 
undetected weak-PAH sources with 2;>1.5. These 
are compared with the three 160 /im -detected weak- 
PAH sources with z> 1.5 (MIPS8342, MIPS15880, and 
MIPS15977). We overlay a number of templates includ- 
ing the classic QSO template of Elvis et al. (1994), a 
torus model, the spectrum of the warm ULIRG Mrk231, 
and an example of a Type-2 AGN template (Lacy et 
al. 2007). We find that the 160 /xm -detected weak-PAH 
sources are roughly consistent with a Mrk231-like SED, 
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Fig. 11. — Top panel: The IRAC8/.tm luminosity vs. Ljr^. 
The solid line is our best-fit relation. We overplot the Bavouzet 
et al. (2008) relations based on their detected sources (red dotted 
line) and including the stacks of the non-detected sources (blue 
dotted line). Whenever Lir is an upper limit, we also give the 
expected range (open circle connected to a filled circle). Middle 
panel: The monochromatic 14 /im luminosity vs. Lj^. The best 
fit linear relation is shown by the solid line. The green dot-dash line 
is the 15 fim conversion relation (Chary & Elbaz 2001). We have 
ignored the difference between the ISO 15 fim luminosity and our 
monochromatic 14 fira luminosity. Bottom panel: The MIPS24 fira 
luminosity vs. Lm- The red dashed and blue dotted lines are as in 
the top panel. In all cases, the red filled symbols are strong-PAH 
sources, while the blue filled symbols are weak-PAH sources. The 
open circles are local ULIRGs (Armus et al. 2007). 



but the stacked 160 /im point of the non-detections (see 
above) is relatively weaker. These are similar to the 
Type-2 quasar from Lacy et al. (2007) (J1711-h5953). 
As is most likely the case with our sources, this source 
is also an AGN-starburst composite (note the presence 
of PAH features). The presence of PAH features after 
stacking the weak-PAH sources (see PaperH), also sug- 
gests that their starburst fractions are > 0. However, we 
should remind that the stacked MIPS160 is only an up- 
per limit. Moreover, half of these sources are undetected 
at 70 /im whose stacked value (1.9 mJy) is about a fac- 
tor of 2 below the faintest of the 70 /im -detected sources. 
Therefore, we cannot exclude that some fraction of our 
sources are consistent with our adopted pure-AGN tem- 
plate. 

Our decomposition implies that the three weak-PAH, 
160 /^m -detected sources have significant starburst com- 
ponents ('--^ 0.5Lir). By analogy, so does Mrk231, which 
has a similar SED. The composite (AGN-|-starburst) na- 
ture of Mrk231 is well estabhshed (Farrah et al. 2003; Ar- 
mus et al. 2007). Based on our AGN-template, Mrk231 
has 63% starburst fraction. By comparison, using the 
Efstathiou & Rowan- Robinson (1995) models, Farrah et 
al. (2003) derive a starburst fraction of 70% for Mrk231. 
More recently, Fritz et al. (2006) derive an total infrared 
starburst fraction of 73% in Mrk231. Our much simpler 
empirical approach comes within 10% of these values. 

To further estimate the range of AGN luminosi- 
ties from different models, we compare our values 



with the SED decomposition of PoUetta et al. (2008), 
which is based on a radiative transfer torus model 
and a starburst template. PoUetta et al. derive 
log Lagn = log ig/im + 0.32. For our sample the median 
(log Lagn - log is.s^iji) is 0.6 which makes Lagn about 
a factor of 2 larger than if we had adopted the Pol- 
letta et al. conversion. However, much of this is due to 
our sample having particularly steep 24 /im /8 /im slopes. 
The PoUetta et al. sample cover a much larger range in 
the 24/im/8/im colors. Their sample, however includes 
5 of our z~ 2 sources (MIPS42, MIPS78, MIPS15840, 
MIPS22204, MIPS22303). Converting their AGN lu- 
minosities to Lq (and accounting for slight differences 
in redshifts), we find that logLACN for these 5 sources 
is respectively: 12.91, 13.07, 12.71, 12.89, and 12.64. 
Our values are always higher, but within <0.1dex of 
the PoUetta et al. luminosities. The only exception is 
MIPS22303 where our Lagn is 0.26 dex higher. These 
examples suggest that our adopted empirical treatment 
of the AGN emission is roughly consistent with other 
approaches; however, the systematic uncertainty of indi- 
vidual sources' Lagn can be as high as 0.3 dex (a factor 
of 2). 

By extension, the 'remainder' of the SED, which is 
attributed to starburst activity is also roughly con- 
sistent. Therefore, we conclude that the starburst- 
fraction of weak-PAH sources ranges from < 30% for 
160 /im -undetected sources (e.g. MIPS78) to ~60% 
(e.g. MIPS15977) comparable to the prototypical warm 
ULIRG Mrk231. 

Accounting for the above systematic uncertainty, even 
the lowest log Lagn of our z > 1.5 weak-PAH sources are 
all still > 10^^ L0. This means that the conclusion that 
our sources host quasar-strength AGN is independent of 
our particular SED fitting approach. 
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Fig. 12. — The average SED of our 160 /im -undetected 
z ~ 2, weak-PAH sources, along with the AGN-component of one 
of our strong-PAH sources (MIPS22482). For comparison, we also 
show the classic Elvis et al. (1994) quasar SED, the Mrk231 spec- 
trum and the SED of a Type 2 quasar from Lacy et al. (2007). All 
data and templates have been normalized at 14 fim . 
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For the quarter of our sample which show strong PAH 
emission features, and typically strong cold dust emis- 
sion (§3.2), the question is reverse from the one in the 
previous section concerning the bulk of our sources which 
are weak-PAH emitters. Based on our spectral decom- 
position, the AGN fraction in the z < 1.5 sources ranges 
from 13% (MIPS22404) to 20% (MIPS283), although 
this could be higher in half the sources due to their 
non-detection at 160 /im. These sources do not show 
signs of either radio or optical AGN. Of the 8 z>1.5 
strong-PAH sources, two (MIPS506 and MIPS15928) 
are actually of borderline PAH-strength (see Paperll). 
Not surprisingly therefore, we find that both show sig- 
nificant AGN fractions (50% and 79% respectively). 
MIPS15928 also shows both a radio and an optical 
AGN (see Table 3). In this section, we concentrate on 
the remaining 6 z>1.5 strong-PAH sources (MIPS289, 
MIPS8493, MIPS16144, MIPS22482, MIPS22530, and 
MIPS22651). Half of these show signs of a radio AGN 
(see Table 3). For these sources our derived AGN frac- 
tions range from 20% (MIPS289) to 32% (MIPS22482). 
Based on our decomposition, the z>1.5 strong-PAH 
sources have Lagn > 10^^ Lq (except MIPS8493, which 
has Lagn = 11.8). These luminosities suggest our sources 
harbor obscured quasars. 

In our fitting, we have assumed that the T ^ 100 K 
continuum (peaking at 20 fim ) is due to an AGN. 
This leads to the 20-30% AGN fractions in our z~2 
strong-PAH sources (see Table 3). However, there is a 
wide range of relative mid-IR continuum levels even in 
pure starburst sources, emphasizing that, especially for 
strong-PAH sources, such infrared SED decomposition 
techniques are uncertain. For comparison, in Figure 13 
we also show the SEDs of the pure starburst, M82, and of 
NGC6240 - a starburst-dominated source but with sig- 
nificant AGN contribution (Risaliti et al. 2006; Armus et 
al. 2007). M82 has the stronger 10-30/im continuum, 
emphasizing the caution needed in assigning this contin- 
uum to an AGN. It is however, worth recalling that M82 
and sources like it are typically much lower luminosity, 
dwarf galaxies - not obvious counterparts of our sources. 

The most closely related population to ours are SMGs 
which have the same redshifts and infrared luminosities 
are our sources. In Figure 13, we compare the spec- 
tra of these 6 z ~ 2 strong-PAH sources with the aver- 
age sub-mm galaxy (Pope et al. 2006, 2008). By av- 
erage sub- mm galaxy we mean, z ~ 2 , 5*850 > 5 mJy and 
Lir~10"Lq. For simplicity, we exclude MIPS22482 
from this figure because its best-fit greybody tempera- 
ture is warmer (60 K) than the rest (~30K). Its AGN 
fraction however is consistent with the three sources 
shown. All our SEDs have relatively stronger < 30 /im 
continua than the SMGs (especially MIPS22482). The 
reality of this can be seen in that even though the red- 
shifts and total infrared luminosities of our sources are 
quite similar to the SMGs, their observed 24 //m fluxes 
(> 0.9 mJy) are consistently higher than those of SMGs 
(Pope et al. 2008). Their observed 70/ini fluxes (^--^5- 
12mJy) are also higher than the stacked 70/Ltm flux of 
z~2 SMGs (-2mJy; Huynh et al. 2007b). Figure 13, 
shows that our strong-PAH sub-population is more sim- 
ilar to local ULIRGs such as NGC6240. This extends 
our analysis of Paperll, where we conclude that the mid- 
IR properties of our sources are very similar to those of 



local ULIRGs, but at higher luminosities. Along with 
the SMGs themselves, this further supports the increas- 
ingly more clear conclusion that the Universe at z ~ 2 
contained starburst-dominated systems of much higher 
luminosities than their local counterparts. 




10 100 
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Fig. 13. — The best- fit SEDs of our strong-PAH 2^2 sources 
(see text for details). These SEDs are compared with the aver- 
age SED of SMGs (Pope et al. 2006, 2008), and the NGC6240 
and M82 templates from the SWIRE template library (courtesy of 
M.PoUetta, see also Weedman et al. 2006). Both our SEDs and 
the comparison templates are normalized so that their 3- 1000 /im 
luminosities are exactly 7 X 10-"^^ Lq . 

3.8.1. Dust mass 

We determine the dust mass only for the sources 
with MAMB01.2mm detections (see Table 1). These 
are most strong-PAH sources (see § 3.8) as well as the 
three weak-PAH sources with strong cold dust compo- 
nent (MIPS8242, MIPS15880, and MIPS22558). We use 
the standard method (see eg. Kovacs et al. 2006), which 
is given in Equation 4. 



(1 + z) X X s(i/rcst,r)' ' ' 

where T is the best-fit cold dust temperature from 
§3.5. The dust absorption efficiency, k is scaled from 
K850 using oc In the literature ksso varies be- 

tween 0.08m^kg~^ and 0.3m^kg~^, depending on the 
dust properties (e.g. size, composition). We use 
Kg50 = 0.15 m^kg"^ as used by Kovacs et al. (2006) for 
a sample of sub-mm galaxies. We obtain dust masses of 
~ 5 X 10* Mq. These dust masses are roughly compara- 
ble to those of SMGs (Kovacs et al. 2006). 

Assuming Mgas/Mdust ^ 100 (Hildebrand 1983), the 
above suggests Mgas ~ 10^" - 10^^ M0. This is compara- 
ble to the gas reservoirs of SMGs. These have been shown 
to be sufficient to sustain the extreme star-formation 
rates needed to power their full infrared luminosisties 
(-lO^^L^) (Tacconi et al. 2006). Initial CO observa- 
tions of some of our mm-bright sources support the spec- 
ulation that these systems also contain large molecular 
gas masses (Tacconi et al. 2008, in prep.). 
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3.9. IR-radio correlation 

The far-infrared luminosity has long been known to 
correlate with the radio luminosity (Helou et al. 1985), 
the ratio of the two usually expressed by q. The original 
definition of q is based on the JKASGO/im and 100 //m 
fluxes, which samples the total 40-120 /xm power (Helou 
et al. 1988). Our data at z^2 spans the restframe 
'--^ 20 - 400 /im , which allows us to measure directly IRAS- 
equivalent luminosity and thus minimize fc-correction ef- 
fects. For consistency, we integrate our best-fit spectra 
to derive L4o_i2o- We rewrite the definition for the q 
parameter to an equivalent, but more convenient form. 

,.,„,(is_H£)_,„,(^),„,3 

To compute the rest-frame 20 cm luminosity, we use 
the spectral indices, a (see § 3.3 and Table 1). For the 
sources without 610 MHz detection, we assume the stan- 
dard a = -0.7. 

The relation remains linear over 4 orders of magnitude 
in luminosity, and for a wide variety of galaxy types with 
g = 2.34 ± 0.01 and a standard deviation of 0.25 dex (Yun 
et al. 2001). The origins of this tight relation are still 
unknown, although likely involve young stars fueling the 
infrared luminosity and, as supernovae, propelling the 
cosmic-ray electrons that fuel the radio luminosity (see 
e.g. Murphy et al. 2006). Early results suggest that the 
relation holds at higher redshifts (Appleton et al. 2004; 
Frayer et al. 2006a). Higher-z can also mean more ex- 
treme luminosities are accessible. Kovacs et al. (2006) 
find that, for a sample of z~2 SMGs, the relation is 
lower than locally (g = 2.07 ± 0.09 using the standard q 
definition) . 
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Fig. 14. — The q-factor (far-IR to radio correlation) as a function 
of redshift. The soUd line represents the local value, and the grey- 
shaded area is the 5cr spread (Yun et al. 2001). As previously, the 
red symbols are strong-PAH sources. The blue symbols and the 
upper hmits are weak-PAH sources. For clarity, we have marked 
the radio-loud sources with open symbols (they are excluded from 
the mean q determination). 

Figure 14 shows the results for our sample. Excluding 
the radio-loud and radio-excess sources (i.e. q< 1.6), we 
find (q) ~ 2.00 ± 0.06 for all sources with radio detections 
(for lower and higher z). There is a weak trend with the 



z ~ 1 (and LiR ~ 10^^ L0) sources being closer to the lo- 
cal relation, while the z ~ 2 (and Ljr, ~ 10^^ Lq) sources 
are consistently too low. In general, steeper radio spec- 
tral indices lower q, while flatter spectra result in higher 
q. To estimate the strength of this effect, we repeat the 
above, however assuming all sources without individual 
measurements of the spectral index have a = -0.5. This 
leads to (g) = 2.04 ± 0.06. Therefore, the uncertainty in 
a does not significantly affect our conclusion of low q val- 
ues. Another potential source of uncertainty is whether 
or not we have underestimated the infrared luminosities 
when fitting our sparsely sampled far-IR SEDs. This is- 
sue was addressed in § 3.6 with the conclusion that this 
is probably not the case and indeed we are more likely 
to have overestimated the infrared luminosities. 

Figure 14 shows that roughly half the sources at z > 1.6 
(excluding the radio-loud sources) are undetected in the 
radio i.e. have only lower limits on q. To estimate the 
effect of this population on the mean q of our sample, 
we use the Kaplan-Meier estimator (Feigelson & Nelson 
1985), which is non-parametric maximum likelihood es- 
timator in the presence of censored data (i.e. lower or 
upper limits). Based on this we revise our estimate of 
for the z> 1.5 sources to (7km = 2.07± 0.01, where (/km 
denotes Kaplan-Meier estimated mean value of q. 

In order to answer the question whether the radio- 
undetected population itself obeys the local relation even 
if the sample as a whole does not, we also stack the 
radio images of the twelve sources without individual 
1.4 GHz detections. We find a mean flux of 57±6.6/iJy 
or an 8 (7 detection. Assuming this value for all twelve 
sources, we find a mean q of 2.21 ±0.02 which although 
still formally lower than the local mean value, is con- 
sistent with it given given the uncertain a values (see 
discussion above). Therefore, the radio undetected 37% 
of the z > 1.5 sources, as well as the z < 1.5 sources are 
consistent with the local IR-radio correlation. 

A lower q than the standard mean value is most com- 
monly but not necessarily associated with excess radio 
emission due to an AGN. However, purely star-forming 
galaxies can exhibit infrared-to-radio ratios a few times 
lower than typical {q^ 2.0) in certain environments such 
as the centers of clusters (Miller & Owen 2001), likely 
as a result of the thermal pressure of the Inter Cluster 
Medium. In rare cases, a merger/interaction can form a 
synchrotron bridge between the merging pair which leads 
to lower q values as well (Condon et al. 2002). When 
considering as yet poorly understood high redshift pop- 
ulations, such effects cannot be ruled out. 

4. DISCUSSION 

4.1. A population of obscured quasars at z ^ 2 

The AGN luminosities we derive based on our spectral 
decomposition ((^agn) — 10^^-^ L0) place our sample in 
the quasar regime. When available, the [Oiii] luminosi- 
ties of our sources, imply intrinsic Ljf-ray ^ 10*^erg/s, 
supporting this conclusion. As discussed in § 3.7 and 
§3.8, our SED decomposition (especially for strong-PAH 
sources) is uncertain. In particular, starburst-powered 
mid-IR continua have been observed both in conjunction 
with strong PAH (e.g. M82) and weak-PAH, but strong 
Si features (e.g. NGC1377; Roussel et al. 2006). The 
presence of either an optical or a radio AGN, however 
is independently supported for at least 20/35 (57%) of 
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our z>1.5 sample (sec Table 3). These include both 
strong-PAH and weak-PAH sources. Although, this docs 
not prove the AGN-origin of their mid-IR continua, it 
does make it the most likely explanation. We therefore, 
conclude that at least 57% of our higher-2; sample are 
consistent with hosting obscured quasars, although this 
fraction can be as high as 100%. The optical faintncss, 
red infrared colors and frequent deep silicate absorption 
all suggests high degrees of obscuration. 

Two sources with silicate absorption do not show opti- 
cal extinction in their Ha/H/3 ratios suggesting perhaps 
an obscured nucleus surrounded by an unobscured host 
galaxy. In Sajina et al. (2007b), we argued that (the 
radio-loud) of these sources might be transition objects 
just before feedback effects dispel the obscuring matter 
and reveal the naked quasar. Some tentative evidence 
of such processes arc seen in the blue-shifted [Oiii] emis- 
sion of MIPS15949 and MIPS16080, which may indicate 
outflows of Au = -576km/s and At; = -261km/s respec- 
tively. 

We find that the number density of our 2; ~ 2 sources 
is ~ 2 X 10~^ Mpc~"^. Their mean AGN luminosity is 
(log(LAGN)) = 12.6 ± 0.2. Since these sources are op- 
tically faint, this luminosity is a close approximation 
to the bolometric luminosity. The number density 
of unobscured (optically and X-ray selected) quasars 
of comparable bolometric luminosities at z^2 is a 
fewx IQ-^Mpc-^ (Hopkins et al. 2007). Therefore our 
sources are already comparable in number to the unob- 
scured quasars of the same luminosity and redshift. The 
mean redshift of our sample, 2; = 2.05, is also consistent 
with the epoch of peak quasar number density which 
is estimated at z^2.15 (Hopkins et al. 2007). Since 
our sample is not a complete census of the unobsciired 
quasars at z ~ 2 (see also Polletta et al. 2008) our results 
support the finding that obscured quasars at ^ ~ 2 are at 
least as numerous, but likely outnumber the unobscured 
ones. 

Lastly, our sample is most likely characterized by 
sources powered by a combination of starburst and 
AGN activity (see §3.7 and §3.8). Concurrent with 
the AGNs arc, frequently ULIRG-strength starbursts 
accounting for ^30-80% of the bolometric emission. 
These support the conclusion that a significant fraction 
of these obscured quasars are accompanied by vigorous 
star-formation activity (see also Martmez-Sansigre et al. 
2008). 

4.2. Implications for the SMG population 

We have shown that 6 of our 35 z > 1.6 sources (20%) 
have 5i.2mm>2mJy, implying S'sso^m ^ 5 mJy. These 
sources therefore would be classified as SMGs in blank- 
sky sub-mm surveys. Their number density however 
(~2/sq.deg.) is only ~ 2% that of the blank-sky sub- 
mm population (Coppin et al. 2006; Pope et al. 2008). 
Four of these six sources are strong-PAH of similar PAH 
strengths as SMGs (Pope et al. 2008). The other 2 
of these strong mm sources, however, have continuum- 
dominated mid-IR spectra (MIPS8242 and AIIPS15880). 
In all cases, the mid-IR continua of our sources are found 
to be a higher fraction of their total infrared luminosi- 
ties than is the case in SMGs (see §3.8). Our model 
suggests that this translates into higher AGN fractions. 
However, as our model assumptions here are uncertain. 



it is more accurate to say that these are upper limits 
on the AGN fractions. Our spectral decomposition sug- 
gests ~ 20 - 30% AGN fractions of Lir for the strong- 
PAH sources, while the weak-PAH sources have 40 - 79% 
AGN fractions. Alexander et al. (2005) find that the 
bulk (75%) of SMGs host AGN; however, even after 
absorption-correction the typical Lx-ray/im, for those 
SMGs is about an order of magnitude lower than that 
of quasars. This implies that their AGN contribute to 
only about 10% of Lir,. This is roughly consistent with 
the estimates of Pope et al. (2007) based on mid-IR 
spectral decomposition. It appears likely that our few 
mm-bright sources represent the stronger-AGN tail of 
the SMG population. Particularly intriguing are the 
mm-bright sources with continuum-dominated spectra, 
because they are true composite objects where AGN and 
star-formation are inferred to each contribute roughly 
50% to LiR. These types of sources are still rare in our ex- 
tragalactic samples, yet are clearly crucial in understand- 
ing any evolutionary links between starburst-dominated 
and AGN-dominated ULIRGs/obscured quasars. Such 
sources are likely to be detected in larger numbers in 
upcoming wide-field sub-mm surveys such as expected 
from SCUBA2 (Holland et al. 2006), when combined 
with large-area mid-IR surveys such as SWIRE (Lons- 
dale et al. 2003). 

5. CONCLUSIONS 

1) Our sample shows low detectability in the far-IR 
(23% at 160 fim, and 16% at 1.2mm). For the ^>1.5 
bulk of this sample, there is a trend of increasing 
MAMB01.2mm fluxes with increasing PAH equivalent 
widths, but no clear trend with increasing 9.7 /um 
silicate absorption. 

2) By contrast 67% of the sample are radio-detected. 
Based on their radio luminosities and/or spectral 
indices, we conclude the 17/35 (48%) of our 2;>1.5 
sources host radio AGN (about half of these are 
radio-loud). This is supported by the observations of 
an enhanced radio emission for a given Lir compared 
with expectations of the local IR-radio correlation. 
Excluding radio-loud sources, and accounting for the 
radio- undetected sources (through survival analysis), we 
estimate (g) = 2.07 ± 0.01, which is significantly below 
the mean value (g = 2.34) estimated locally. This radio 
excess is most likely, but not necessarily, due to the 
presence of radio AGN. 

3) We further test for the presence of AGN using 
optical/near-IR spectroscopy which is available of 21 of 

the sources. We find a good correspondence between 
the mid-IR and optical diagnostics in that only 10% 
of the strong-PAH sources show signs of optical AGN 
versus 62% of the weak-PAH sources with optical/near- 
IR spectra. This supports on average the use of 
PAH-strength as a discriminant of AGN-dominated 
vs. starburst-dominated sources. However, 5 of our 
weak-PAH sources show starbm-st-like optical spectra. 
Two of these show signs of radio AGN. 

4) Of the 11 sources with near-IR spectra, we find 
only one broad line source 

(FWHM(Ha) >2000km/s). This is in contrast with 
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Brand et al. (2007) who find that 7/10 of their high-z ob- 
scured ULIRGs show broad emission lines. On the other 
hand, the optical spectra of four of our sources show the 
Civ emission line with widths corresponding to the NLR. 

5) Evidence for blueshifted [Oiii] emission is seen in 
two of our z ^2 wcak-PAH sources. This is potentially 
indicative of outflows of up to 570km/s. 

6) Where available, we find no correspondence be- 
tween optical extinction and the depth of the silicate 

absorption feature. We estimate that the lower limit 
on the star-formation rates in z ^ 2 weak- PAH sources 
without far-IR detections is 2O-3OM0/yr. By contrast, 
far-IR detected, z^2 sources (typically strong- PAH, see 
below) have SFR~ QOOMo/yr. 

7) We fit the available far-IR data with a simple grey- 
body model. We combine this with our mid-IR fitting, 
in order to estimate the total starburst Imninosity (cold 
dust plus PAH), and AGN luminosity (hot/ warm dust 
continuum) . In sources with high PAH equivalent widths 
roughly ^70 80% of their total infrared luminosity is 
due to starburst activity. In weak-PAH sources this is 
typically < 30%. More intriguing arc a a few weak-PAH 
sources with far-IR detections, suggesting significant 
starburst contributions (-^50%). 

8) We confirm expectations, that our sources have 
bolometric infrared luminosities of ~ 10^^ Lq for the 
2; ~ 1 sub-sample, and close to ~ 10^'^ Lq for the 2; ~ 2 
majority of our sample. We present relations between 

and i24Aim and Lir. 

9) Compared with SMGs, our strong-PAH 2~2 
sources typically have similar redshifts, total lumi- 
nosities, sub- mm emission (and hence dust mass), but 
higher mid-IR continua (attributed to greater AGN 
fractions). Our SED fitting suggests these sources 
have ~20-30% AGN contribution to their Lm- The 
presence of AGN is supported in 3/6 of these sources 
which have AGN- like radio luminosities. For our z^2, 
strong-PAH, MAMBO-detected sources, we find dust 
masses of ^ 5 x 10* M0, consistent with the values for 
sub-mm galaxies. This suggests that these sources 
have comparable gas reservoirs to SMGs and hence are 
capable of sustaining the SFR^ 900 M0 /yr implied by 
their AGN-subtracted infrared luminosities. 
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Far-IR-to-radio properties of Spiteer-selected high-^ ULIRGs 



TABLE 1 

FAR-IR and radio fluxes with 1 (T UNCERTAINTIES. 



MIPSID 


z 




5l60'' 


Si. 2mm 


S20cm 


•SeiOMHz 


_,610MHz 
"l.4GHz 






mjy 


mJy 


mJy 


mJy 


mJy 


strong-PAH 


283 


0.938= 


6.4±2.2(2.8") 






0.16± 0.02 


0.35±0.09 


-0.96 


289 


1.86 


5.6 ±2.1(3.6")'' 


54±10(20.0")'' 


2.28±0.35 








506 


2.470^^ 






1.37± 0.53 


0.14± 0.03 




>-0.92 


8184 


0.99 


17.9±1.8(3.0") 


49±9(5.0")'' 


0.72±0.62 


0.16± 0.02 


0.23±0.05 


-0.43 


8207 


0.834« 


17.7±1.5(0.7") 


48±8(4.8") 




0.28± 0.03 


0.50±0.06 


-0.72 


8493 


1.80 






0.94±0.52 




0.18±0.05 


<-0.83 


15928 


1.505<= 


12.4±1.9{1.2") 




-0.04±0.56 


1.25± 0.06 


2.93±0.15 


-1.02 


16144 


2.132<= 






2.93±0.59 


0.12± 0.03 




>-1.10 


22404 


0.610= 


47.1±4.8{3.1")':: 


134±12(9.4")'' 


0.66± 0.59 


0.20± 0.03 




>-0.49 


22482 


1.84 


12.0 ±1.5(0.1") 


34±7(1.5") 


1.46±0.59 


0.25± 0.02 


0.61±0.09 


-1.05 


22530 


1.952': 




38±9(9.0") 


2.11±0.56 


0.15± 0.03 


0.49±0.07 


-1.45 


22554 


0.824': 


10.2 ±1.7(2.2") 






0.13± 0.03 




>-1.01 


22600 


0.850^: 


6.4±1. 7(6.9") 












22651 


1.73 


5.5 ±1.2(3.8") 


37±12(9.8") 


2.18±0.54 


0.47± 0.03 


0.94±0.09 


-0.83 


weak-PAH 


42 


1.95 


10. 7±1. 3(3. 0) 




0.43±0.d2 








78 


2.65 






-0.25± 0.62 








110 


1.050^: 


9.5±2.2(7.2")'^ 




-0.42±0.55 


0.37± 0.03 


0.81±0.16 


-0.94 


133 


0.905=: 


12.3±2.0(1.3") 




0.01±0.57 


0.39± 0.03 


0.73±0.09 


-0.77 


180 


2.47 


5.8±1. 7(4.3") 




1.26±0.44 








227 


1.635': 


10.1±2.0(3.4") 




-0.41±0.58 


0.30± 0.03 




> -0.00 


279 


0.952<: 


9.5±2.1(1.0")'' 


27±9(11.7")'' 


0.31± 0.50 


0.35± 0.06 


0.86±0.19 


-1.10 


429 


2.213=: 






1.03±0.57 








464 


1.85 


7.3±1. 5(3.2") 




-0.27±0.65 


0.15± 0.02 




>-0.86 


8034 


0.95 






0.74±0.62 


0.14± 0.02 




>-0.92 


8196 


2.586=: 


5.0±1. 4(5.6")'' 




0.99±0.43 








8242 


2.45 


4.8±1. 2(1.3") 




2.52±0.52 








8245 


2.70 






-0.49±0.46 


0.130± 0.009 


0.24±0.06 


-0.74 


8268 


0.80 






-0.45±0.78 


0.047± 0.008 






8327 


2.445=: 






1.03±0.59 


1.40± 0.06 


3.44±0.14 


-1.08 


8342 


1.562=: 


10.7±1. 4(0.9") 


29±11(8.8")'' 


0.98±0.52 


0.18± 0.03 


0.86±0.11 


-1.90 


15840 


2.30 






-0.08±0.57 








15880 


1.68 


19.0±2. 1(1.1") 


23±8(9.3") 


2.06±0.53 


0.60± 0.21 


1.35±0.06 


-0.98 


15949 


2.118': 


7.3±1.6(0.6") 




1.24±0.51 


0.16± 0.02 




>-0.76 


15958 


1.97 






1.28±0.50 


0.27± 0.02 


0.39±0.07 


-0.45 


15977 


1.85 


14.3±1.0(1.1") 


40±7(3.5") 


1.38±0.53 


0.298± 0.008 


0.61±0.07 


-0.95 


16030 


0.98 






0.09±0.69 








16059 


2.325=: 






1.20±0.66 


0.57± 0.03 




> 0.77=1 


16080 


2.007=: 


5.2±1. 7(6.0") 




0.69±0.54 


0.34± 0.03 


0.63±0.07 


-0.75 


16095 


1.68 


11.0±1.6(1.3") 




-0.65±0.50 








16113 


1.90 






0.72±0.52 








16122 


1.97 






0.00±0.60 


2.82±0.14 


4.34±0.12 


-0.52 


22204 


1.974=: 


12.9±1.8(2.7") 




0.27±0.51 


1.89±0.08 


4.67±0.10 


-1.09 


22277 


1.77 


17.5±1.8(2.2") 




0.86±0.6 


1.57± 0.07 


4.79±0.11 


-1.34 


22303 


2.34 


6.6±1. 5(1.4") 




-1.05±0.64 


0.26± 0.02 




>-0.17 


22467 


0.70 


9.1 ±1.5(3.3") 




-0.73±0.70 








22558 


3.20 






1.73±0.48 


0.62± 0.03 


1.55±0.12 


-1.11 


22661 


1.80 






0.65±0.56 








22699 


2.59 






-0.31±0.63 




0.25±0.07 


<-1.23 



^We use a format of Flux it Icr. The offset between the APEX source and 24;Ltm source is given in brackets.. Large 
offsets are typically blended sources (see 
'^Multiple 24 jLtm sources. We give the deblended value (see §2.10.1 for details). 

^Redshift from optical or near-IR spectroscopy (see §3.4). These are estimate to be accurate to ±0.001. The rest of 
the redshifts are based on mid-IR features (see Paper I). 
Potentially an inverted spectrum (see §3.3 for details). 
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TABLE 2 
3(T Detection statistics 




All 


strong-PAH 


weak-PAH 






z<1.5 


2> 1.5 


z< 1.5 


z>1.5 


70 fim 
160 fim 
1.2mm 
1.4 GHz 
610 MHz 


60% (29/48) 
23% (11/48) 
16% (7/44) 
67% (32/48) 
48% (23/48) 


100% (6/6) 
50% (3/6) 
0% (0/2) 
83% (5/6) 
50% (3/6) 


50% (4/8) 
50% (4/8) 
50% (4/8) 
75% (6/8) 
50% (4/8) 


57% (4/7) 
14% (1/7) 
0% (0/7) 
71% (5/7) 
43% (3/7) 


56% (15/27) 
11% (3/27) 
11% (3/27) 
59% (16/27) 
48% (13/27) 
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TABLE 3 

Fit results 



MIPSID 


log(Ll.4GHz) 


log(LAGN)'' 


log(LsB)'' 


log(LlR)<='d 


T'cold 


SFR 


Radio 


Optical 




W/Hz 


Lq 


Lq 


Lq 


K 


M@/yr 


AGN 


AGN 


strong-PAH 


283 


23.84it0.16 


11.3±0.1 


<11.9 


11.3-12.0 


40±5 


< 110 


no 


no 


289 


<24.2 


12.1±0.3 


12.7±0.2 


12.8±0.2 


29±2 


880±220 


no 




506 


24.67±0.14 


12.7±0.2 


<12.7 


12.7-13.0 


34±5 


< 730 


no 


no 


8184 


23.75±0.23 


11.6±0.1 


12.3±0.2 


12.4±0.2 


43±2 


320± 90 


no 




8207 


23.89±0.08 


11.4±0.1 


12.1±0.2 


12.2±0.2 


40±1 


210± 60 


no 


no 


8493 


<24.2 


11.8±0.3 


<12.4 


11.8-12.5 


32±4 


< 450 


no 




15928 


25.27±0.01 


12.6±0.1 


<12.1 


12.6-12.7 


25±4 


< 200 


yes 


yes 


16144 


24.46±0.43 


12.3±0.2 


<12.8 


12.3-12.9 


30±3 


< 920 


no 




22404 


23.43±0.16 


ll.litO.l 


11.9±0.2 


12.0±0.2 


35±2 


130± 40 


no 


no 


22482 


24.81±0.06 


12.4±0.2 


12.8±0.2 


12.9±0.2 


60±3 


900±360 


bor 




22530 


24.83±0.22 


12.4±0.3 


12.8±0.2 


12.9±0.2 


33±2 


990±400 


bor 


no 


22554 


23.54±0.29 


11.0±0.2 


<11.8 


11.0-11.9 


44±5 


< 100 


no 


no 


22600 


<23.4 


11.1±0.2 


<11.7 


11.1-11.8 


39±4 


< 90 


no 


no 


22651 


24.92±0.03 


12.2±0.1 


12.7±0.2 


12.8±0.2 


30±2 


860±210 


bor 




weak-PAH 


42 


<24.3 


13.0±0.1 


<12.4 


13.0-13.1 


30±5 


< 380 


no 


— 


78 


<24.6 


13.1±0.1 


<12.7 


13.1-13.2 


36±6 


< 730 


no 




110 


24.33±0.04 


12.0±0.1 


<11.9 


12.0-12.2 


23±3 


< 110 


no 


yes 


133 


24.13±0.04 


11.7±0.1 


<11.8 


11.7-12.1 


43±5 


< 90 


no 


yes 


180 


<24.5 


12.9±0.2 


<12.6 


12.9-13.1 


35±5 


< 700 


no 




227 


24.61±0.06 


12.6±0.1 


<12.2 


12.6-12.8 


28±4 


< 290 


yes 


yes 


279 


24.24±0.08 


11.6±0.1 


11.9±0.2 


12.1±0.2 


40±3 


140± 60 


no 


no 


429 


<24.4 


12.7±0.2 


<11.7 


12.7-12.8 


20±0 


< 90 


no 


no 


464 


24.42±0.40 


12.5±0.2 


<12.1 


12.5-12.7 


26±7 


< 220 


no 




8034 


23.72±0.09 


11.9±0.1 


<11.8 


11.9-12.1 


21±2 


< 90 


no 




8196 


<24.5 


12.9±0.1 


<12.7 


12.9-13.1 


39±5 


< 880 


no 


yes 


8242 


<24.5 


12.7±0.2 


<12.9 


12.7-13.1 


34±3 


<1160 


no 




8245 


24.76±0.05 


12.8±0.1 


<12.7 


12.8-13.1 


39±6 


< 790 


no 




8268 


<23.4 


11.4±0.1 


<11.6 


11.4-11.8 


20±0 


< 70 


no 




8327 


25.88±0.02 


12.7±0.1 


<12.7 


12.7-13.0 


36±5 


< 750 


yes 


no 


8342 


24.82±0.08 


12.4±0.1 


12.1±0.2 


12.6±0.2 


27±2 


220± 70 


bor 


yes 


15840 


<24.4 


12.7±0.1 


<12.7 


12.7-13.0 


37±4 


< 770 


no 




15880 


25.06d=0.03 


12.8±0.1 


11.9±0.1 


12.9±0.1 


20±0 


140± 30 


yes 




15949 


24.60d=0.20 


12.8±0.1 


<12.5 


12.8-13.0 


33±4 


< 560 


no 


yes 


15958 


24.62±0.11 


12.5±0.1 


<12.5 


12.5-12.8 


32±3 


< 520 


no 




15977 


24.81±0.08 


12.6±0.1 


12.8d=0.2 


13.0±0.2 


56±2 


920±230 


bor 




16030 


<23.6 


11.5±0.1 


<11.8 


11.5-12.0 


22±2 


< 110 


no 


no 


16059 


25.23±0.04 


12.7±0.1 


<12.6 


12.7-12.9 


33±5 


< 610 


yes 


yes 


16080 


24.88±0.05 


12.6±0.1 


<12.3 


12.6-12.8 


30±6 


< 320 


bor 


no 


16095 


<24.1 


12.6±0.1 


<12.2 


12.6-12.7 


27±5 


< 230 


no 




16113 


<24.2 


12.2±0.2 


<12.5 


12.2-12.7 


32±4 


< 470 


no 




16122 


25.68±0.01 


12.4±0.2 


<12.5 


12.4-12.8 


32±4 


< 480 


yes 




22204 


25.78±0.01 


12.9it0.1 


<12.4 


12.9-13.0 


31±5 


< 380 


yes 


yes 


22277 


25.69±0.01 


12.8±0.2 


<11.8 


12.8-12.9 


20±0 


< 90 


yes 




22303 


24.90±0.09 


12.9±0.1 


<12.5 


12.9-13.0 


32±6 


< 500 


yes 




22467 


<23.2 


ll.litO.l 


<11.5 


11.1-11.7 


39±4 


< 50 


no 




22558 


25.83±0.02 


13.3±0.2 


<12.0 


13.3-13.3 


20±0 


< 160 


yes 




22661 


<24.2 


12.2±0.1 


<12.4 


12.2-12.6 


31±4 


< 420 


no 




22699 


<24.5 


12.6±0.1 


<12.8 


12.6-13.0 


39±4 


<1060 


no 





^The integrated 3 — 1000 fim luminosity of the AGN component alone. 

^The AGN-subtracted Lib.- Includes both cold dust and PAH emission. It is an upper limit for all sources without 

160 i-tm. detections. In these cases, the 2 a limits arc used in the fits. 
*^The total integrated 3— 1000/im luminosity. 

■^For 160 //m -undetected sources, we take Lib. to be in the range from Lagn (i-e. mid-IR continuum) alone to 
^AGN and the upper limit on LsB- 
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TABLE 4 
Near-IR spectroscopy results 



MIPSID zopt FWHM(Ha) log(LHc) FWHM(Oiii) log(Loiii) log(Nii/Ha) log(Oiii/H/3) Classification 







km/s 


Lq 


km/s 


Lq 








429 


2.213 


918 


8.61 






<-1.00 




Starburst 


506 


2.470 


743 


8.38 






-0.46 




Starburst 


8196 


2.586 


562 


8.84 


702 


9.02 


-0.91 


0.69 


Composite 


8327 


2.441 


971 


8.92 


634 


8.94 




0.54 


Composite 


8342 


1.562 


771 


7.61 






-0.18 




Composite 


15949 


2.122 


2189 


9.21 


1304 


8.93 


-1.44 


> 1.06 


AGN 


16059 


2.326 


666 


8.21 


671 


8.56 


-0.64 


> 1.70 


AGN 


16080 


2.007 


1359 


9.13 


1155 


8.88 


-0.88 




Starbursf 


16144 


2.131 


522 


8.68 






-0.40 




Starburst 


22204 


1.974 














Blended?'' 


22530 


1.952 


461 


8.37 






-0.86 




Starburst 



^The line ratio is starburst-like. However, this source shows the asymmetric [OIII] profile typical of AGN (see §3.4.3). 
'^This source is classified as a Type 2 AGN based on its UV spectrum (Martmez-Sansigre et al. 2008). 



